Quantum irreversible decoherence behaviour in open quantum systems with few 
degrees of freedom. Apphcation to NMR reversion experiments in nematic hquid 

crystals. 



H. H. Segnorile and R. C. Zamar 

Instituto de Ftsica Enrique Gaviola - CONICET, Facultad de Matemdtica, 

Astronomia y Ftsica, Universidad Nacional de Cordoba 
M.Allende y H. de la Torre - Cmdad Umversitaria, X5016LAE - Cordoba, Argentina 

(Dated: May 7, 2013) 

An experimental study of NMR spin decoherence in nematic liquid crystals is presented. The 
outcome of the experiments are analyzed in the framework of a theory that considers the spins as 
an open quantum system coupled to a quantum molecular environment, presented by the authors 
recently. Decoherence dynamics can be put in evidence by means of refocusing experiments of 
the dipolar interactions. The experimental technique used in this work is based on the MREV8 
pulse sequence. Non-idealities of the experimental setting, like external field inhomogeneity, pulse 
misadjustments and the presence of non-reverted spin interaction terms are analysed in detail and 
their effects on the observed signal decay are estimated. It is found that, though all these non- 
idealities could in principle affect the evolution of the spin dynamics, their influence can be mitigated 
and they do not present the characteristic behaviour of the irreversible spin decoherence. As unique 
characteristic of decoherence, the experimental results clearly show the occurrence of eigen-selectivity 
in a time scale intermediate between those controlled by quantum interference of the closed spin 
system and thermalization, in complete agreement with the theoretical predictions. We conclude 
that the eigen-selection effect is the fingerprint of decoherence associated with a quantum open spin 
system in liquid crystals. Besides, these features of the results account for the quasi-equilibrium 
states of the spin system, which were observed previously in these mesophases, and lead to conclude 
that the quasi-equilibrium is a definite stage of the spin dynamics during its evolution towards 
equilibrium. 



I. INTRODUCTION 

Derivation of a thermodynamic-like stationary state, 
of equilibrium or quasi-equilibrium. from a mieroseopic 
quantum mechanical starting point, is a major prob- 
lem faced today by the physics of irreversible processes 
in nonequiUbrium systems, with impact in a variety of 
fields, from fundamental research to applications in ar- 
eas of high current interest like quantum computing and 
quantum information theory 

Particularly, the occurrence of quasi-equilibrium spin 
states in nematic liquid crystals (LC's) poses the problem 
of the irreversible evolution of a finite open quantum sys- 
tem of interacting particles, coupled with a large quan- 
tum environment, towards a quasi-stationary statc^^. In 
LC's, due to the rapid molecular motions that average the 
intermolecular spin interactions to zero, the effective spin 
system comprises a small number of magnetic degrees of 
freedom, namely the intramolecular dipolar interactions, 
which remain because of the typical oricntational order of 
these mesophases^^. However, the proton Nuclear Mag- 
netic Resonance (NMR) response in LC's is consistent 
with a true quasi-equilibrium in spite of the small num- 
ber of the degrees of freedom of the observed system^'^. 
This topic has attracted the interest of many researchers 
in the NMR field "'^^""'^^ as well as other areas of physics 

In our recent publication on NMR quantum decoher- 
ence in LC's, which we shall refer to as QD-I^^, a the- 
ory was presented which describes the irreversible quan- 
tum decoherence processes undergone by an observed and 



controlled system of quantum interacting particles be- 
cause of its coupling with an unobserved lattice or envi- 
ronment. It provided a comprehensive explanation, com- 
patible with previous experimental evidence^^'^*^, on the 
mechanisms which turn the proton spin system density 
matrix of nematic LC's into a quasi-equilibrium form af- 
ter an arbitrary initial coherent state . 

Accordingly, the decoherence process transforms the 
initial state into a block diagonal matrix in the cigenbasis 
of the spin-environment interaction Hamiltonian. This 
evolution occurs over a time scale intermediate between 
that of the Liouvillian evolution of an isolated observed 
system and the long time scale where evolution is gov- 
erned by relaxation and thermalization processes driven 
by thermal fluctuations. 

The fingerprint of quantum decoherence theory is the 
eigen-selection process, which causes a selective decay of 
the off-diagonal components of the density matrix, pre- 
serving the block diagonal part of the initial state. Such 
explanation relies on general requirements on the eigen- 
value distribution functions of the relevant quantum op- 
erators that represent the interaction of the system with 
the environment. Hence, experimental observation of 
eigen-selectivity would provide a direct evidence of the 
correlated dynamics between the spin system and the 
quantum environment that preludes the reach to quasi- 
equilibrium. 

The aim of this work is the experimental study of deco- 
herence spin dynamics, within the theoretical framework 
of QD-I. Particularly, we focus on the most relevant as- 
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pects of this phenomenon, namely eigen-sclcctivity and 
the occurrence of the intermediate time scale. Irre- 
versible decoherence is studied by means of refocusing 
experiments, designed to counteract the effect of the spin 
dynamics generated by the dipolar spin interactions. Be- 
sides of the quantum interference corresponding to a 
closed system, the intrinsic decoherence coexists with 
other sources of the signal decay, like the line broaden- 
ing due to a distribution of the order paramcter^"'^, non- 
idealities of the experimental setting like inhomogeneity 
of the static and rf mangetic fields or pulse imperfec- 
tions, and the effects of non-secular terms of the dipolar 
Hamiltonian that cannot be experimentally rcversed^^. 
Accordingly, to isolate the distinctive features of decoher- 
ence associated only with the microscopic dynamics, we 
explore the intermediate time scale by combining dipolar 
refocusing with a meticulous theoretical and numerical 
analysis of the experiments based on hypotheses of gen- 
eral character. This strategy enables us to visualize the 
physical processes involved and to quantify the influence 
of the non-idealities. 

Confirmation of the irreversibility of the spin dynamics 
within the intermediate time scale has relevant practical 
and basic consequences. It implies that the state attained 
by the spin system over the intermediate scale has a rep- 
resentation in the form of a block diagonal density opera- 
tor, containing the information of the initial preparation 
which does not undergone subsequent evolution, except 
for the slow dynamics imposed by the spin-relaxation 
process, no matter the details of the coherences present 
in the initial condition. The traditional introduction 
of the spin-temperature hypothesis amounts to assum- 
ing that the off-diagonal density matrix elements can be 
forgotten after quantum interference cancels their con- 
tribution to the observable expectation value^^. In this 
work we provide a quantum explanation for the damp- 
ing of the off-diagonal elements and for the time scale 
where it occurs. Clarifying this aspect enables to replace 
the phenomenological assumption by a specific condition 
satisfied by the spin system state. Also, understanding 
the nature of the decoherence mechanism can contribute 
to the set up of a quantum spin-lattice relaxation the- 
ory beyond the Markovian limit, providing insight into 
the interplay of the quantum correlations developed in 
the microscopic irreversible dynamics and the dissipative 
macroscopic evolution. Finally, the analysis carried out 
in this work may contribute to the current discussion 
about the role of system-environment entanglement as 
basic mechanism of quantum decoherence of interacting 
particle systems^'^^"^^. 

The article is organized as follows. In Section II, by 
means of experiments of refocusing of the dipolar inter- 
actions, we show how the eigen-selection process is evi- 
denced during the partial reversion of the spin dynam- 
ics, which occurs in an intermediate time scale, between 
the time scales of the Liouvillian free evolution or adi- 
abatic by one hand, and the thermalization and relax- 
ation processes by the other hand. The initial part of 



Section II is devoted to a detailed summary of the con- 
cepts developed and results obtained in QD-I, which are 
exhaustively used in this work. In Section II A a theoret- 
ical description of the decoherence dynamics driven by 
reversion experiments is presented, while a series of ex- 
perimental measurement are shown in Section II B and 
the results compared with the theoretical prediction. Sec- 
tion II C is concerned with the analysis of possible sources 
of eigen-selectivity like effects coming from a closed-spin- 
system dynamics under experimental misadjustments or 
non-idealities in the theoretical approach used for its de- 
scription. A supplementary matcrial^^ is attached to this 
work, where the details about the origin of the theoret- 
ical expression for the reversion evolution operator and 
the signals obtained in Section II B are brought together 
the results of numerical calculations for such signals un- 
der different experimental settings. Besides, an experi- 
mental analysis of the effects of the inhomogeneity of the 
static magnetic field is presented there. 
Finally, Section III is devoted to discuss and conclude 
about the concepts developed and results obtained along 
this work. 



II. EIGEN-SELECTIVITY EFFECTS ON 
COHERENCE EVOLUTION 

This section is devoted to the experimental detection of 
the eigen-selection process which characterizes the irre- 
versible time evolution of the spin coherence in LC's, due 
to the coupling of the spin system to a quantum lattice 
or environment. This will be done within the framework 
of the theoretical proposal presented in QD-I. 

For the convenience of the reader, hereinafter we in- 
clude a brief summary of the steps followed in QD-I to 
work out the nonequilibrium irreversible dynamics of an 
observed and controlled system coupled to an external 
unobserved environment. This problem was approached 
using the following Hamiltonian: 

H = Hs + Hsl + Hl, (1) 

where Tis = 1^^^ and "Hl = l'"^ n^/^ are 

the Hamiltonians of the observed system and the lat- 
tice or unobserved environment respectively. Symbols of 
the form O*^") and O^-'^) indicate operators acting exclu- 
sively on the Hilbert space of the system and the lattice 
respectively. The interaction Hamiltonian acts on both 
Hilbert spaces, and is represented by Hsl = J2q^q-^q^ 

with Fq = ® F^-^\ and A, = A^"^ <Si Index q can 
label different characteristics, like a spin pair or a tensor 
component of the interaction Hamiltonian. 

The existence of distinct time scales associated with 
different physical processes, which become important as 
the spin dynamics evolves, was postulated. The occur- 
rence of such time scales is subjected to certain general 
conditions that the relevant Hamiltonians must fulfill. 
Such rcqiiircments arc associated with commutation re- 
lations between T-Ls, T-Ll and Hsl (it is worth to note 
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that [HstHl] = 0). In this way, the dynamics of the 
observed system in the earher time scale is obtained by 
assuming that ['Hs,'Hsl] = and [Hlt'Hsl] — 0, which 
impUes a system model called essentially isolated system, 
where the observed system evolves reversibly and the ex- 
pectation values of the observables decay due to quantum 
interference. A later time scale is obtained by assum- 
ing [Hst'Hsl] = and [?^l,?^sl] 0, where the sys- 
tem model was referred as essentially adiabatic system, 
and the corresponding dynamics of the observed degrees 
of freedom is irreversible. Finally, the latest time scale is 
obtained by assuming ['Hs,'Hsl] and ['Hlj'Hsl] ^ 0, 
and the system model was named system in thermal con- 
tact. 

The focus in QD-I was laid on the study of pure de- 
coherence processes of the observables due to the cou- 
pling of the system to a quantum environment, with- 
out thermalization and relaxation effects. Therefore, the 
dynamics of the observed system was studied assuming 
[HsjUsl] = 0, namely, the system dynamics was de- 
scribed under essentially isolated or essentially adiabatic 
models, where the observed system conserves its energy. 
Within this framework, NMR experiments on a wide va- 
riety of nematic LC's can be described using the following 
Hamiltonian definitions: 

Us = Y,Usi, HsL = YlHsLi, (2) 

i i 

where 

■Hsi = -Wo Izi + ("Hdi)/, (3) 

H5L. = H(f ®(si{)-5,,l(/)). (4) 

Index i labels the i-th molecule (the sum runs over the 
molecules of the whole sample), uq = jBq is the Larmor 
frequency, 7 the proton gyromagnctic ratio and Bq the 
strength of the static magnetic field, which is applied 
along the laboratory z axis, and l^i is the z projection of 
the spin angular moment for all the protons in the i-ih 
molecule. 

Following, we will describe the nature of the Hamilto- 

nians in Eq. (2) . The operator H^^) in Eq. (4) is the spin 
part of the contribution Us Li to the interaction Hamil- 
tonian HsL between the observed system (the protons) 
and the environment, and is defined as 




The indices j, k run over all the proton sites within the 
z-th molecule, rjk is the internuclear distance between 

spins 7 and fc, fijk stand for the polar angle of the vec- 
tor rjk with respect to the system fixed to the molecule. 



The spin angular moment is Ij = Ixj x + lyj y + Izj Z- 
The introduction of the quantum character of the envi- 
ronment variables through the operators S^^i ^1- (^) 
is a novel contribution of QD-I, which allows a general 
description of the dipolar Hamiltonian in liquid crystal 
NMR with the important consequence that the quantum 
correlation between the observed system and the envi- 
ronment can be included in the dynamic. Precisely, as 
shown in QD-I, the irreversibility of the spin dynamics 
naturally emerges in the decoherence time scale when 
the full-quantum character of the spin-environment in- 
teraction energy is assumed. The meaning of S^^^i is 
that of an operator whose eigenvalues S^zi are the possi- 
ble values of the zz component of the order parameter, 
Szzi = (1 c^os^ ^2^5); where 6i is the angle between the 
long molecular axis of the i-th molecule and the external 
static magnetic field Bq. Then, the full-quantum dipolar 
Hamiltonian of nematic LC's is written as 

nd = Y,Udi, (6) 

i 

where 

^., = HW®si{i, (7) 

with 

H|jf = l^^^i) (g) • • • (g) H|f'^ (g) • • • 

where the superscript (si) indicates an operator acting on 
the spin Hilbert space of the z-th molecule. In the def- 
inition of Eq. (6), a motionally averaged dipolar Hamil- 
tonian is used, in which the intermolecular spin interac- 
tions are averaged out due to the fast molecular diffusion. 
This means, in a quantum language, that the eigenval- 
ues of the intermolecular lattice operators of the dipolar 
Hamiltonian arc negligible in comparison with the in- 
tramolecular ones and the time scale of the dynamics 
originated from them is much longer than that of the 
intramolecular spin interaction. Accordingly, such slow 
dynamics can be neglected. The average dipolar energy 
is obtained through the lattice density operator as 

{Hdi)f = trj {Hdi PL{eq) } , (8) 

where the trace is taken over the lattice variables, (•)/ = 
trj {•PL(eg)} is the expectation value, and PL{eq) is the 
lattice density operator at thermal equilibrium, 

PLie,) = l^'^^e-^-<"/Aff, (9) 

where /3t = fe^) the Boltzmann constant, T the 

absolute temperature, and A// = trf ^e~^'^^'''^^^. The 
dipolar Hamiltonian corresponding to the closed spin sys- 
tem is obtained after the trace over the lattice variables 
is taken in Eq. (8): 

= 5,, H(f®l(/), (10) 
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with S,, = (S[{h 



Finally, using the Zeeman Hamiltonian 

Hz = — WqIz = —Wo y^Izi; 



(11) 



and the average dipolar Hamiltonian of the spin as a 
closed system 



(12) 



we define the Hamiltonians in Eq. (2) as 

Hs = Hz + {Hd)f, (13) 

and 



nsL = nd-{Ud)f. 



(14) 



In closing the definitions of the Hamiltonians of 
Eq. (1), it is worth to mention that the environment 
or lattice Hamiltonian 'Hl is associated with the po- 
tential energy of the mechanical interaction between the 
molecules of the whole sample. This Hamiltonian takes 
into account the molecules as an ensemble of correlated 
quantum objcts. The details of the mechanical inter- 
molecular interaction arc not needed to extract some 
general conclusions about the influence of the environ- 
ment on the spin dynamics. In fact, knowing whether 
T-Ll commutes or not with the interaction Hamiltonian 
HsL enough for concluding about the reversible or ir- 
reversible character of decoherence in some time scale, 
as was shown in QD-I. However, the influence of Hl is 
indirectly taking into account when wc introduced the 
molecular averaged values of lattice operators, as is seen 
in Eq. (8), and over the distribution probability function 
of the eigenvalues of such operators, as we will show in 
Section II A. 

In the following, we will discuss how decoherence in- 
duced by the environment emerges under the dynamics of 
the system described by the Hamiltonian of Eq. (1). We 
are concerned with expectations values of spin observ- 
ables, therefore, we need to calculate the reduced spin 
density matrix. First, we introduce the eigenbasis of the 
operators H.^^) , {| C-s ) = | Ci^i ) ® • • • (S) | Qsi ) (g) • • • (g) 
I CivSjv)}, which span the spin Hilbert space of the A'" 
molecules of the sample. The symbols Q^s indicate the 
different eigenvalues and Si label their degeneration, thus 

■^di'l = Cj| C*)- Then, we calculate the density ma- 
trix time evolution under the Hamiltonian of Eq. (1). It 
is convenient for comparison with NMR experiments to 
use the rotating-frame representation, namely a frame 
whose z-axis is parallel to the external magnetic field, 
and rotates about it with an angular frequency equal to 
the Larmor frequency. Thus, in this frame the time evo- 
lution of the density matrix is 



p(t) = U(i)&(0)pL(e«)UHi) 



(15) 



where O = e ^'^"i^* O e*""'-^* is the representation of op- 
erator O in the rotating-frame and 



Vit) 



(16) 



where Hs = J2i{'^d.i)f is the transformed spin Hamil- 
tonian. The initial state of the spin system ps{0) is ob- 
tained, for example, after applying a radiofrcquency pulse 
sequence. Therefore, the matrix elements of the reduced 
density operator, a{t), in the basis {| )} are 



{Cs\d{t)\Cs')^trf{{Cs\m\Cs')} 



(17) 



where S ^ (H^fV 



Szz and 



%,C'}W = trf {ut(/)(C',i)U(/)(C,t)pi^(i,J , (18) 

is the decoherence function associated with the time evo- 
lution operator 



U(/)(C,t) = e- 



(19) 



In order to explicit the spin dynamics of a representative 
molecule, say the i-th one, we use the operator expansion 
technique due to Zassenhaus to factorize the evolution 
operator (19) into a product of exponential operators, as 
followsii'28 : 

U(/)(C,t) = e-*««'(';)*e-*';'«^'"*ug'^^^^^(C,f), (20) 
where we defined the operator 



■n 



if) 



(21) 



In turn, the third factor in (20) can be expanded as 



U 



imsLi) 

C 



(22) 



where C 



imsLi) 



{() represents some ((;-l)-order nested 



commutator between Hamiltonians v}Q{C,) and CiU'^sLi- 
In keeping with the same spirit of separating the molec- 
ular spin dynamics, we consider the cases where the ini- 
tial condition /Os(0) can be expressed as 



&(0) 



1 

W- 



(23) 



where Nsi = trs^ {l^'^^^} is the trace of the identity op- 
erator in the Hilbert space of the spins belonging to a 
molecule, pij^^O) = l'^^) «) • • • ® p("')(0) ® ■■■ ® 
where p*^*'^(0) is the initial condition of the spin vari- 
ables of the i-th molecule, which is assumed identical for 
all molecules. Then, the matrix elements of (23) in the 
spin space are 

( wm cs' ) = E ( \psim > 

' (24) 



=^E(c.«^i^"^(o)ic^«on^§:S' 
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where we defined 

as the product of the TV — 1 Kroncckcr deltas associated 
with all the molecules different from the i-th one. Then, 
by using Eqs. (20) and (24), Eq. (17) can be expressed 
as 

{Cs \a{t)\ Cs' ) = J2e-'^^^-^'^^''{Qs, |p (-)(0)| ({s^ ) 

i 

(25) 

where G{f^^'}(i) is the decoherence function (18) for the 
cases where the Kronecker deltas present in Eq. (24) do 
not cancel. Because of the delta functions, (18) is evalu- 
ated when C,j = Vj /j^i for the expression shown in 

Eq. (20), hence n^^}{C) = n^nHO and in lj[^'>^'^^'\c,t) 
only Q is unaffected by such delta functions (i.e. it can be 
Q ^ Q). This reasoning leads to the following expression 
for the decoherence function 

^ ' . (26) 

xug)(-^)(C,.)pi^(i,,}, 

where we introduced the operator \j''^^^f''"'\c,,t) which 

can differ from Up ((^, i) only by the replacement of 
Ci by Ci- It is worth to remark that if Q = (I in (26), 
then the decoherence function satisfies 

G{U.}it)=trf {Pi?e,)} = l' 

which indicates that decoherence does not affect the sub- 
space associated with a given degenerate eigenvalue Q. 
An important consequence of this is that, if in the ma- 
trix representation of the spin density operator p^^*\0) 

in the eigcnbasis of T-l^gl^ the eigenstates are ordered, in 
such way that blocks associated with a given eigenvalue 

are formed in the diagonal, these blocks remain invari- 
ant under decoherence. This property, which was already 
reflected in Eq. (18) for ^ = C is of great importance in 
the dynamics that brings the spin system into a quasi- 
equilibrium state. 

Finally, a zero-trace spin observable acting on the spin 
spaces of individual molecules, and whose operator is 
time independent in the rotating-frame, has the general 
form: 

6 = ^6f)®l(/), (27) 

i 

with of-' = 1(^*1) (g) • • • (g) O^**') (g) • • • (g) By using 



Eq. (25), the expectation value of O can be written as 

(28) 

with 

G{C.,caW = '^^E^ ' G{,,,,y{t)/Ms^-\ (29) 

where the sum in Eq. (29) runs over all the values of 
and Sj. with k ^ i. 

The irreversibility of the spin dynamics is introduced 
through of decoherence associated with an open quan- 
tum system characterized as an essentially adiabatic sys- 
tem. According to the theory, the eigen-selection process, 
prominently involved in the decoherence process, leaves 
invariant the block diagonal space of the density operator 
drawing it to a quasi-equilibrium representation. 

As will be shown below, the experimental measure- 
ments exhibit the eigen-selection effects during the co- 
herence evolution under a process of partial reversion of 
the spin dynamics. Besides, two different time scales as- 
sociated to different sources of decoherence will be distin- 
guished, one of them during free evolution and the other 
under time reversion of the coherences. In Section II A 
the reversion pulse sequences used in the experiments are 
described. The experiments arc designed for the case of 
evolution of single quantum coherence corresponding to 
the FID signal. However, the procedure and the results 
can be extended to any kind of coherence. The measure- 
ments and the analysis of the results are presented in 
Section II B. Finally, in Section II C the different sources 
of error in the experiments are considered. 

A. Definition of the experiments 

In order to highlight the contribution to the coher- 
ence signal decay coming from sources other than the 
Liouvillian dynamics corresponding to an isolated spin 
system (i.e. quantum interference), we will measure the 
single coherence signal evolution of the FID, when the 
spins are subjected to a rf pulse sequence configured to 
compensate the Liouvillian evolution. An experiment of 
reversion of the spin dynamics under the high field sec- 
ular dipolar Hamiltonian can be performed by means of 
the sequence known as MREV8^^"^^, which is shown in 
Figure 1 (b). An explanation of the effect on the spin 
dynamics due to such kind of reversion sequence can be 
found in the supplementary material (see supplementary 
materiaP^, section "Misadjustments in the dynamics un- 
der reversion" , for details on the reversion spin dynam- 
ics) . This sequence was selected because that the quality 
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in the reversion process rehes on the relationship between 
the total time of the FID evolution and the setting time 
between pulses. In this reversion pulse sequence, the 
smaller the setting times n and T2 are, the better the 
reversion will be. Therefore, to access to long reversion 
periods while keeping the time parameters as small as 
possible, we apply a pulse sequence consisting of a chain 
of blocks of the same reversion unit, as is shown in Fig- 
ure 1 (c). On the other hand, in other techniques like 
the 'magic sandwich'^^'"^^"'^^ the mentioned quality relies 
on the high intensity of the rf pulses and it could become 
inadequate for cases of long time evolutions due to the 
requirement of rf pulses of long duration, as is the case 
of the nematics VKAdQ and PAA studied in this work. 

For the sake of simplicity in the analytical calculation, 
we will first use the simple unit of Figure 1 (a) to ob- 
tain the analytical decoherence function under reversion 
of a FID. This procedure will allows us to extract general 
conclusions which can be later extended to the sequences 
used in the experiments, such as those shown in Figure 1 
(b). To study decoherence under reversion, the spin dy- 
namics can be considered as a composition of a free evo- 
lution during ti followed by a 'reversion' evolution due 
to the two pulses separated by the time t^, as shown in 
Figure 1 (a). By applying the free evolution operator 
(16) over the state | Qs), it is obtained 



(30) 



where U(-^)(C,<) was defined in Eq. (20). Truncating the 
expansion of Eq. (22) to the first order gives 



(31) 



which is valid for the intermediate time scale, under 

the postulate of the existence of different time scales 
described in the introduction of Section II. The anti- 
Hermitian operators 

cSl(C) ^ [EO^si'^s'^J, (32) 



p(/) 



ray(/) -iy(/) 1 



(33) 



have pure imaginary eigenvalues. 



The evolution operator of the dynamics under rever- 
sion in the rotating frame can be written, in general, as 
(see section III.D.2 of QD-I): 



Urt(t) 



i Hs t/K i (Ei nsLi/K -Ul) t 



(34) 



where k is a positive constant whose value depends on 
the particular refocusing technique used. The evolution 

operator factor g*^"'*/'' for the closed spin system un- 
der an MREV8 experimental reversion sequence, with a 



particular time setting, is obtained in the supplementary 
material (sec supplementary materiaP^, section "Misad- 
justments in the dynamics under reversion" , for details 
about the analytical expression of the evolution operator 
of the MREV8 reversion dynamics). For the cases (a) 
and (b) of Figure 1, k = 2. By applying (34) over the 
state I Qs ) and truncating as in Eq. (31), we get 

Urt (0 1 Cs ) = ^' I ) ® U(^') (-C/« , t) , (35) 

with 



(36) 



which is also valid for the intermediate time scale. 
i f) 

n^-{—(/K) is obtained by replacing Q by —Q/n, Vj, 
in Eq. (21). 

The complete evolution operator of the reversion se- 
quence of Figure 1 (a) will be the composition of Eqs. (16) 
and (34), that is 



U(t2,il)=Urt(i2)U(il), 



(37) 



written in the rotating system. For a single sequence 
we have ti = ti and t2 = ^2, with ri and T2 defined in 
Figure 1 (a) . In the case of the chaining block sequence 
shown in Figure 1 (c), the complete evolution operator 
has the same form as Eq. (37) except that ti = nri and 
t2 = nT2, being n the number of blocks. Again, by ap- 
plying (37) to the state | )) we obtain 

U(i2,ii)| Cs) = e-*^>^>^(*-*^/-)| Cs) ® U(^nC,i2,ti), 

(38) 

where 

U(/)(C,i2,ii) ^ e-^^'^-^/-^)*^ e-''^«''(^)*i 

X e-'i''^sL e'^M^sLiO (t?-Ft^/«"-2tit2/«)/2 

eC*c<f2 (t?-t^/«-2tit2/«)/2 



X 



(39) 



Eq. (39) represents the dynamics under an MREV8-like 

sequence when the experimental setting has ri and T2 
small enough to neglect the non-secular terms of the 
dipolar Hamiltonian which arises from applying the 7r/2 
pulse pairs"^^ (see also supplementary matcriaP^, section 
"Misadjustments in the dynamics under reversion", 
for details about the influence of such non-secular 
dipolar terms on the MREV8 reversion dynamics). This 
constraint is imposed over each block in the sequence of 
Figure 1 (c) but the total times nri and nT2 do not need 
to be so small. 

Finally, the evolution operator that represents the spin 
dynamics due to the pulse sequence of Figure 1 (c) is 
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completed by adding in Eq. (37) the time operator cor- 
responding to the FID evolution during the time t, that 
is 



U(t,i2,ii) = U(i)Urt(t2)U(ii), 
with the corresponding effect over the state | (s ) 



u{t,t2,h)\Cs) = 



— p-^Y.^C^S{t + t^-t2/K.)\ 



where wc defined 

U(/)(C,t,t2,ii) ^ U(^''(C,i)U(^HC,i2,ti), 
which can be approximated by 

U(^) (C, i, t2 , i 1 ) = e-^ ' *e-* *^ 



(40) 



(41) 



(42) 



(43) 



X e 



Ci^ifi (f^+2tt2+2tti+ti-t2/K-2tif2/K;)/2 



for small ti and r2 as in Eq. (39). Now, using Eq. (43) 
we can obtain, with the same spirit as in QD-I, the de- 
coherence function for the spin dynamics valid for a first 
time scale longer than the Liouville characteristic time 
scale, that is 

xU(/)(C,t,i2,ii)pi^(i,)} 

KCi-C^)iC<^]j- (C) {t^-2tt2/K+2ttl+tl+tl/K.^-2tlt2/K)/2 



X e 



X e 



(Ci-C,')iC<f2 {f^+2tt2+2tti+tl-tl/K-2tit2/K)/2 (/) 



(44) 

where we introduced the operator 'UiP{(^,t,t2,ti) which 
can differ from U^-'''(C,t, t2,ii) only by the replacement 
of Ci by Ci- The operators iC^^^ and iCyi appearing in 
Eq. (44), are Hermitian due to the anti-Hermitian char- 
acter of the commutators, as we have seen in (32) and 
(33), thus their eigenvalues are pure real numbers. By 
making t2 = Kti (i.e. T2 = kti), we eliminate the Liou- 
ville dynamic during the reversion and Eq. (44) becomes 
in this case 



X e-'(fi-^.')^<^'f2 {[*+(«+!) *il'-(«'+3«+2)t?}/2 (45) 
X p-i(f*-0'C<fijC)tV2 „(/) \ 



with G{i^^Q]{t,ti) = G{(;_(;'}(f, Kti,ti). Notice that 
Eq. (45) depends just on the reversion time ti, due to 
the election t2 = Kt^. 

To extract a final expression for the decoherence func- 
tion, we use that the environment states form a continu- 
ous and dense space, so that we can replace in the trace 
the sum by an integral over the lattice space, obtaining 



^-Ca^*'"^) =/ j J '^^ ^f'9K,h P^/ec,) 
X e-'(Ci-C:)ASi(/i)tg-i(c,-c.')cf^''('»)tV2 (46) 
X e"' iii-OcHa) {(t-Fr)=-[l-F(«-M)-i] r^}/2^ 

where r = {k + \)ti is the total time under rever- 
sion, kf^g = (/Is), kg^h = {9\h) and p^'^^^^ = 
{^\('^L(eq)\ ^ ) ■ have used in (46) the eigenbasis 

{|/)}, {Is)} and of the Hermitian operators 

U^il, iC\2 andiC,%(C), respectively, withW^^^J /) = 

\.fi) ® ■■■ ^n'ii^lfi) ■■■ ^ I/n) = A^,(/,)| /), 

A5,(/,) = S,if,)-S,,, ^Ci2\9) = C[ig)\g) and 
*clgi(C)| /i> = C^^'^{h)\h), where ASi{fi), C[{g) and 
C^^{h) are real numbers. Also, we used the closure 

relations Jdg\g){g\ = l'-^^ and Jdh\h){h\ = I'^^l 
Now, defining the eigenbasis of the environment or lat- 
tice Hamiltonian Hl as {|/^)}, with the eigenvalues 

Ej{f^), where | ) = Ej (f^) | ) and us- 

ing the closure relation J df^\ )(f^ \ = we can 
write: 



Pl'U = I df' kf.j pj^^,,) (/^) , (47) 



where 



PL{eq) 



The final decoherence expression is obtained changing 
in Eq. (46) the state integration variables df , dg and dh, 
by the eigenvalue integration variables rfAS",, dCf' and 
dCf^, thus we can finally write 



X J dC^ e-*(«'-«>)^f {(*+^)'-[i+('=+i)"']^'}/2 (48) 
X J dCf,^ e-'^i^-Ocft*'/^ p,iASuC^,Cfl), 
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where it is defined the function 

X j df j dg j dh kj^g kg^t k^jE kfEj 



dr 

Ef,AS, (49) 



and [■]^L^sL indicates that the integrals over d/^, df , 

dg and dh are calculated maintaining constant the eigen- 
values Ej\ ASi, C[ and Cf^, respectively. Since 

Pi{ASi,Cl,Cfl-) satisfies 

J dAS, J dC^ J dCfl'pi{ASi,C^,Cf^) = 1, 

it can be interpretated as a probability distribution 
function of the eigenvalues given by ASi, Cf' and Cf^. 

We can see from Eq. (48) that the decoherence function 
is the result of a superposition of complex exponential 
functions weighted by a distribution of their frequencies. 
Such kind of superposition can be generically expressed 
as follows: 



Vi{ACi,t,T) = dAie 



I' 



Pi{Ai,ACi,t,T), 
(50) 

with AQ = Q — and yj^. (t, t) is some polynomial in 
t and T, the distribution of frequencies of the complex 
exponential is determined by Pi{Ai, A(i,t,T). The ex- 
pression (50) is a general form for the different integrals 
in Eq. (48). Therefore, we can use (50) to extract con- 
clusions about (48). 

We are concerned with studying the decoherence function 
produced by the coupling of the spin system with an en- 
vironment whose states belong to a continuous and dense 
Hilbert space, where a distribution of the eigenvalues of 
each complex exponential in Eq. (48) can be defined. In 
the physic systems we are concerned with, such distribu- 
tion is supposed to have a general bcU-shapc form around 
the mean value of the eigenvalues, converging to zero fast 
enough so that integrations in Eq. (48) and Eq. (50) can 
be extended to ±oo. Prcssumably, these conditions on 
the eigenvalue distribution functions are similar to those 
imposed in the work of reference^^ to derive the equili- 
bration of 'quasi-isolated quantum systems' in the strict 
sense. 

We can observe that Eq. (50) is the Fourier transform 
on the variable Ai of the function Pi{Ai, AQ, t, r) valued 
in A(i {t, r), that is 



Vi{AC,i,t,T) =TAi {Pi(^i,ACi,i,T)}^^.3;^^(j_^) 



/+(X> 
-oo 



dAi e 



-i OL Ai 



Pi(yli, ACj,t,r) 



Therefore, Eq. (48) is the Fourier transform over all the 
variables of pi(AS'i, C/", C^^), which is valuated in some 
polynomial in t and/or r multiplied by A<Ji for each vari- 
able. Due to the general bell-shape form assumed for 
the distribution functions Pi{ASi,C^ ,Cf^) characteriz- 
ing the environment of our physical system, the super- 
position of complex exponentials functions which consti- 
tutes (48) and (50) will have a form of a decay function 
in the time t and/or r. It is easy to see that the bigger 
the value of A^j is, the faster the decay of such func- 
tion will be. In the case of A^j = (i.e. = C) we 
have G'{(^.(^'}(f, r) = 1 and we do not have a decay. All 
these characteristics are described as the eigen-selectivity 
or the eigen-selection effect over the dynamic of the ob- 
served system, as we have seen in QD-I and at the end 
of the introduction of Section II. 

This behaviour of the decoherence function, produced 
by the eigen-selectivity, is still obtained from a non- 
truncated time evolution operator, having a more com- 
plex form than (43). In such case the form of the de- 
coherence function will be like (48) but with a more ex- 
tensive development of integrals due to a bigger amount 
(maybe infinite) of eingenvalues coming from nested com- 
mutators of increasing order, and with a more complex 
distribution function pj. Therefore, a complete form (i.e. 
without trrmcating) of the decoherence function (48) can 
be expressed as 



'''''' D,n 



dAS^ e-^^^-'^^-'pi ( AS.. 



(52) 



with 



(51) 



where JJJ Ylj:, „{'} represents multiple integrals over 
the eigenvalues of different class of nested commutators 
labeled by D and with a nesting order n. Such eigen- 
values are generically termed C/^'"". The polynomials 

yD.n(t,T) are functions of t and/or r which multiply 
the corresponding Cf/" eingenvalue in the complex 

exponential functions. The symbol |c/^'"| repre- 
sents a set of such eingenvalues. It is worth to note 
that the general expression (52) is valid for any kind 
of experimental setting, the only difference between 
experiments can be the forms of the polynomials 3^n,„ 
with a possibly different time dependence. 

To extract a more handleable version for the decoher- 
ence function, we will use the hypothesis of the existence 
of different time scales in the dynamics, as was assumed 
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previously in the introduction of Section II. The suitabil- 
ity of this assumption to describe the physical systems of 
interest was experimentally confirmed in other works of 
the authors^^'^'^. Under such hypothesis, the dynamics 
produced by the exponential operator with n'gli, which 
was called reversible adiabatic quantum decoherence in 
QD-I, is faster than the dynamics produced by the re- 
maining exponential operators with nested commutators, 
which was called irreversible adiabatic quantum decoher- 
ence in QD-I. Hence, the decay produced by the com- 
plex exponentials with ASi is faster than the ones pro- 
duced by the remainder exponentials with eigenvalues 
^D,n rpj^g experimental evidences^ ^'^^ show that these 
two time scales are very well separated. Therefore, in 
Eq. (52), for a time t over which the decay function pro- 
duced by the complex exponential with A^^ is close to 
vanish, the values of the exponentials with eigenvalues 
Cj^'" do not significantly deviate from their values for 
t = 0, thus the decoherence function can be very well 
approximated replacing the multiple integrals over such 
eigenvalues by their values in t = 0, namely, it can be 
approximated using /// JJj^ ^ {-j^^^. 
Finally, applying this approximation to the particular 
case of Eq. (48) gives 



with 



-i(C.-C')AS«t 



G{CcMt^r)c^ J dAS, 



(53) 



where it is defined the function 



= j dEf PL(eq){Ef) 

j df^ j df j dg kf^g kgje kfEj 



(54) 



Ef,ASi 



From Eq. (53), we see that under such approximation the 
dependence of the decoherence function on Cf^ is elim- 
inated in Eq. (48), thus the decoherence function only 
depends on the values of Q and Q. Therefore, we will 
have from Eq. (29) 

Two additional hypotheses of general character will be 
assumed on the decoherence function, as follows: 

HypoG I: Statistical independence between the 
eigenvalues ASi and C/", Vi. This hypothesis is 
reasonable due to the different nature of the spec- 
tral properties of their operators. This is written 

as Pi{ASi,Ci^) - ^> (ASi) pf"^ (Cf) and it 
brings us the possibility of writing Eq. (53) as 



= ydA5,e-'(^*-O^^**pP>(A50, (55) 

and 

g|J^,j(t) ^ J dCf e'^^^-Oo^r-/[^i-+^)] pf''} {Cf) . 

(56) 

HypoG II: Homogeneous environment for each 
molecule, i.e. absence of border effects. This 
is written as Pi{ASi, Cf ) = pj{ASj,Cf), and 
it brings us the relationship 

At this point, some comments about the MREV8 
pulse sequence used in the experiments are pertinent. 
This sequence is composed of two sequences WHH-4^^, 
which correspond to each half in Figure 1 (b). Due to 
the symmetric disposition of the pulses in the AIREV8 
sequence, it can be seen that the dynamics under 
reversion does not depend on the first-order nested 



(C) in Eq. (22), which means that 

(/) 



commutator Q^^^^^^^^ 

the dynamics under reversion produced by C^l in the 
decoherence function (45) is completely reverted when 
h = K^i- Therefore, in such reversion experiment, 
the decoherence function with the evolution operator 
truncated up to the first-order commutator will be equal 
to the decoherence hmction (45) putting ti — 0, and it 
will be equal to (46) or (48) putting r = 0, so they will 
be independent of the reversion time. 
Obtaining the decoherence function under MREV8 
using a truncated expression for Eq. (22) up to the 



second-order conmutator C 



(C) would involve a 



G 



tough calculation. 

However, the approximations and conclusions extracted 
from the general decoherence function form (52) will still 
be valid and the result obtained in Eq. (53) will differ 
in the complex exponential function involving Cj^. This 
function will have a dependence with other additional 
eigenvalues C^^'" and their corresponding polynomials 
yD,n{t,T). Since that aim is out of the scope of our 
work, we will conserve the decoherence function (53) 
for our succeeding analysis. Due to the hypotheses 
HypoG-I and HypoG-II will be valid in the general 
case and the final form of decoherence depends on a 
proposed distribution function of some eigenvalues (as 
we will sec the in the following text), the conclusion will 
be unchanged if we use the asymmetric single reversion 
sequence of pulses shown in Figure 1 (a). 

We will mention too that the small ti approximation 
can be improved by adjusting the time between pulses 
and the pulse width as small as possible. As was men- 
tioned, the factor k in the operator Eq. (34) is equal two, 
thus ideally the sequence MREV8 of Figure 1 (b) will 
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revert the spin dynamics when the condition T2 = 2ti 
is satisfied. The reversion experiment for the FID with 
MREV8 wiU be that shown in Figure 1 (c) wliorc each 
block RS(n) is a sequence Uke the presented in Figure 1 
(b). This configuration allows to set the time between 
pulses as small as allowed by the experimental appara- 
tus, in order to minimize the effects of non-idealities of 
the sequence. In this case, the time period of evolution 
under time reversion is t = riTc, where n is the number 
of MREV8-blocks, and the FID is acquired during time 
t. 

The observed FID corresponds to the expectation value 

of the observable ly, ensuing the reversion sequence. 
From the theory presented in QD-I, which was revised 
in the beginning of Section II, using Eq. (28) with 



p(«i)(0) = ^4'^^ and 6(^1) 



r(si) 



and the deco- 



herence function (53) under HypoG-I and HypoG-II, 
in the 'on-resonance' condition, we have: 



Ci^ijCi^i 



(57) 



= c7 = /«(/) \ = 



where 5 = = (S, 



zzl 



= Szz and the values of r are 



multiples of Tc- 

The result of Eq. (57) is similar to the FID obtained in 
Section III.E in QD-I, with the addition that the signal 
is attenuated by a decoherent factor given by g'"^*^^, (t) , 

which depends on eigenvalues Ci and C,[ of the dipolar 
Hamiltonian. 



In Eq. (57), the function G(^_^_(^'^{t) corresponds to the 
main decoherence process under a free evolution of the 
system. This decoherence was called adiabatic quan- 
tum decoherence (AQD) in QD-I, with a time scale 
shorter than the rest of the decoherence processes and 
a reversible character of its dynamics. As was ana- 
lyzed in QD-I, in nematic liquid crystals the AQD is 
associated with the orientational molecular distribution 
function (OMDF). That distribution is expressed as a 
distribution of the values of the order parameter Si. 
The distribution function associated here with such de- 
coherence is p\^^^ i^Si), being ASi = Si — S and 
S = J dSi Si p\^^ (Si). Using the random variable trans- 
formation (RVT) theorem^^, we have that p\^^ (Si) = 
^{AS} _ ^j^g distribution function pj^^ is the 

same that the one given by pj^^^ but with its statisti- 
cal variable shifted in the mean value S), thus the AQD 
function in (57) can be expressed as 



/oo 
-00 

/OO 
dS, e-*(^^-«i)^^*pP (5i). 

(58) 

From Eq. (58), we can see that G{(;^^j}(t) is the same as 
the one obtained in Section III.C in QD-I, which exposes 
the relation between the AQD function and the OMDF 
given by p{'^^ (^i). 

In QD-I, we have seen that a gaussian distribution is 
a suitable approximation for the OMDF in nematic, this 
is 



1 



2ir cr|^ 



Si 



(59) 



with asi the standard deviation of Si which is the same 
for A^i. 

Finally, wc can sec that the distribution (59) yields a 
gaussian form for the AQD function (58) namely 



(60) 



In Eq. (57) , the function Gj.^*^, (r) corresponds to the 
decoherence function produced by the first-order nested 
commutator belonging to a set of nested commutators 
which conform the complete evolution operator. The dy- 
namics produced just by this decoherence function could 
be eventually reverted, by using a particular pulse set- 
ting, but it is impossible to revert simultaneously the dy- 
namics of the whole set of commutators, as was demon- 
strated in Section III.D.2 of QD-I. Therefore, such deco- 
herence function produced by the complete set of nested 
commutator is irreversible; it was called essentially adi- 
abatic quantum decoherence in QD-I and it has a slower 
dynamics than the reversible AQD, introducing an in- 
termediate time scale between the AQD or Liouvillian 
process and the process of thermalization. 

We can obtain different decay functions for G^^*^, (r) 

depending on the form of p\^ ^ i^i) ■ ^'^^ instance, if 
we have for such distribution a gaussian form 



(61) 



thus decoherence has the following gaussian behavior 

<)(r) = e-^^-^^'^^^r^/[«(«+^)1. (62) 



By other hand, for a Lorentzian distribution 



(63) 
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decoherence has the following exponential decay behavior 



(64) 



In the following, we will calculate the Fourier trans- 
form on the time t of Eq. (57) with the aim of analyzing 
the spectral properties of the former function. First, we 
know that the Fourier transform of a complex exponen- 
tial function is a Dirac delta function, thus we obtain 

{e-*(^i-Ci)^*} (a;) = 27r5 [cv - {([ - C,)S] . (65) 

Next, we will calculate the Fourier transform of the AQD 
function, that is 

/C50 
-OO 

/oo 
dAS, e-^«i-«i)^^^*pi^^> (A5i) 
-OO 

/oo 
dASi S [CO - id - Ci)ASi] p{^^> (ASi) , 
-CSO 



(66) 



where we have used that 
S 



[co - id - Ci)A5i] = ^ r dt e-' [— (Cl-Ci)ASi] 

Using in Eq. (66) the delta function property 

1 



5[a;-(Ci-Ci)A5i] 



ICi-Cil 



(5[A5i-c./(C;-Ci)], 



we have the following result for the Fourier transform of 
the AQD function 

The last step to calculate the Fourier transform in the 
time t of Eq. (57) is to obtain the Fourier transform of the 
product of a complex exponential function with the AQD 

function, whi(;li c;an be; calculated using the convolution"^^ 
of the Fourier transforms (65) and (67), thus we obtain 

j-^{e-(Ci-c;)s*Gc^_c;(i)}(a;) 

= {e-K-^i)^*} (c.) * M 



X'l -_Ci 

_ 27r {AS} ^ c.-(C(-Ci)^ 



Pi 



(68) 



where A{ijj) * B{lS) is the convolution between the func- 
tions A and B, and we used in (68) that 



Finally, using (68) in (57), we obtain the Fourier trans- 
form of the signal produced by the reversion experiment 
as 



J-*{(ly(i,r))}(a;) 



fi«i.f(«i 



(69) 



2^ plAS} ^-jCl-Cl)^ ^) ^(.*)^(^)_ 



la-Ci 



a-Ci 



At this point, we will analyze the meaning of the result 
in Eq. (69). First, we note that if G{^^^^/}(f) = 1 and 

^c''c' ^'^) "^(^j ''')' (^''')' ^'^^^ obtain the resulting 
signal of a FID under a description of the spin system as 
a closed system. Therefore, the Fourier transform of such 
signal will be a superposition of Dirac deltas shifted at 
the spectrum characteristic frequencies (Ci — C,i)S, like 
(65). By other hand, if in (57) and (69) we just assumed 
G^^*^,{t) ~ 1, Vr, wc would obtain the resulting signal 
of a FID under the AQD dynamics, which is the same 
result obtained in Section III.E of QD-I. 

The AQD produces the line-shape of the spectrum. 

Such spectrum is obtained as a superposition of copies of 
the OMDF shifted at the frequencies (Ci— Ci)"^ scaled 
by the factor \([ — Cil) ^ can be seen in (69). This scal- 
ing effect of \<^[ — Cil produces the cigen-selection effect 
in the time domain, due to which the bigger the value 
of ICi ~ Cil is, the faster the decay of the decoherence 
function will be. Such effect is reflected in the spectrum 
making that for higher frequencies {([ — (i)S the corre- 
sponding copy of the OMDF is less intense and wider. 

Now we gained more insight to understand the result- 
ing signal under reversion in Eqs. (57) and (69). We can 
see that the spectrum (69) resembles the one obtained 
for the AQD in in QD-I, but in the present case there is 
a modulation produced by decoherence during the rever- 
sion dynamics. Such decoherence under reversion pro- 
duces a faster decay if the value of \([ — Cil is higher, as 
can be seen from Eqs. (62) and (64). Therefore, if we 
have different spectra for different values of r we should 
see that the higher the frequency (Ci — Ci)*? of the spec- 
trum line is, the faster the decay in r will be, provoking a 
kind of compression of the spectrum. This feature of the 
dynamics under reversion will provide us a clear method 
to detect the effects of the eigen-selectivity due to the 
coupling of the system with the environment. 

As a final comment, from Eq. (51) and Eq. (56), we can 
observe that the decoherence function under reversion is 
the Fourier transform in the variable Cf" of the distribu- 

tion function p] , valuated in —{Q — ^•) r /[2(/t -|- 1), 
this is 
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f 

.J — c 



2(^ + 1) 



' 2(« + l) 



(70) 



The expressions (67) and (70) show the close relation 

existing between the dccohcrencc process and the 
distribution function of the variables associated with the 
environment. 

The numerical calculation of Eq. (57) is shown in Fig- 
ures 2 and 4 for 5CB and PAA^g respectively, where a 
frequency selective gaussian decay was proposed for deco- 
herence along t and r time scales. In the case of 5CB, we 
used the ten-spin model (core protons plus the a — CH2 
proton pair) from reference^^ because of the high effort 
that would involve the inclusion of more spins, which 
is beyond of the current computational facilities. For 
PAAd6, the eight-spin model presented in QD-I was used. 
In Figure 2 (a) it can be seen the detail of the FID signals 
for different values of t in 5CB, while Figure 2 (b) shows 
the corresponding amplitude spectra. The frontal view 
of the spectra in Figure 2 (c) shows that the faster de- 
cays correspond to the components of higher frequency. 
This feature is more evident for the normalized spectra 
of Figure 2 (d), where it can be appreciated the com- 
pression of the spectrum for increasing values of r, in 
correspondence with the smoothing of the signals as seen 
in Figure 2 (a). 

The details of the variation of the normalized ampli- 
tude of the spectrum lines of frequencies 5.65 KHZ and 
8.50 KHz is shown in Figure 3, where it can be appreci- 
ated the higher decay rate of the high freqiiency compo- 
nent. The calculations presented in Figure 2 show the ef- 
fects on the dynamics introduced by the eigen-selectivity 
of the decoherence process. The choice of the 5CB sample 
for demonstrating the effects of the eigen-selectivity was 
motivated by the clear separation existing between the 
spectral peaks of the two groups of spectral lines around 
5.65KHZ and 8.50KHz found in the numerical calcula- 
tion, corresponding to the strong dipolar couplings of the 
molecular core and the a-pair respectively. On the con- 
trary, the spectrum of PAA^g does not exhibit a clear 
distinction among the frequency lines, since the strong 
dipolar couplings have similar values. However, a com- 
pression of the spectra as a function of the reversal time 
T similar to that found in 5CB is still visible in this com- 
pound, as can be seen in Figures 4 (a) and (b), because 
of which it is used to show that this effect is part of the 
evidence of the eigen-selectivity of the decoherence pro- 
cess. Notice that if decoherence did not present eigen- 
selectivity during the reversion period, all the spectral 
lines would decay with the same rate and therefore the 
spectral compression would not occur. 

Finally, it is worth to remark that the occurrence of 



decohercnt factors is a consequence of incorporating the 
'mechanical' variables into the description of the dy- 
namics within a full-quantum view, and that the eigen- 
selectivity is a direct consequence of this fact. As will be 
shown in Section II C, the relevant experimental errors 
do not entail eigen-selectivity. The experimental mea- 
surements corresponding to the experiments proposed 
are presented in Section II B . 



B. Measurements 

In this section we present the experimental results ob- 
tained by application of the refocalization experiment of 
Figure 1 (c) in the FID dynamics, which was discussed 
in detail in Section HA. 

The experiments were carried out in a home-built spec- 
trometer, based on a magnet of a Varian EM360, of 60 
MHz for protons, with the probe adapted for application 
of pulsed radiofrequency. The electronic setup allows the 
complete control of the phase of the pulses, with a pre- 
cision of 0.022°, a time step of 40 ns, and a minimum 
configurable time of 240 ns. The rf power used permits 
rf 7r/2 pulses of about 6.5/iS. The homogeinity of the 
magnetic field is controlled with shimming coils, with 
a minimun half width of 470Hz approximately, for the 
spectral lines (using a model of gaussian line-form for 
each spectral line). The maximum dead time in the sig- 
nal acquisition is about 18/is. The temperature can be 
set between 25°C y 150°C, with a medium accuracy of 
il^'C and a stability of ±0.1°C. 

We made experiments in samples of the ncmatic liquid 
crystals 5CB (4'-pentyl-4-biphenyl-carbonitrile), PAA^g 
(methyl deuterated para-azoxyanisole) y PAA (para- 
azoxyanisole) , and the solid adamantane. The solid sam- 
ple was included for comparison, since this system con- 
tains an unlimited number of interacting spins, in con- 
trast with the finite 'clusters' comprising the protons of 
LC molecules. Degradation of spin coherence could also 
occiir in a solid induced by non-spin dc^grcxs of freedom, 
for instance due to lattice phonons. However, this point 
is out of scope of the present work. 

With the purpose of reference for subsequent discus- 
sion, in Figure 5 we show the measured FID's with the 
corresponding spectra, for all the samples studied. These 
signals are the result of eight acquisitions, except in 
PAArfg, where the signal was acquired 208 times. PAA^jg 
and PAA molecules differ in that in the first the methyl 
groups are replaced by CD3 groups. In adamantane, 
while the molecules as a whole are fixed on the solid 
network, they undergo rapid motions, due to which the 
crystalline spin system can be effectively represented by 
a lattice of spins 1/2*'^. 

It should be noted that the experimental spectrum of 
5CB in Figure 5 (a2) is different from the spectrum nu- 
merically calculated with the ten-spin model used in Sec- 
tion II A. It can be seen that in the measured spectrum 
the higher amplitude peaks are shifted to higher frequen- 
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cies, while in the calculated spectrum the situation is the 
opposite. This is so because the ten-spin calculation does 
not include the remaining CH2 proton pairs of the alkyl 
chain, which contribute to the high frequency peak. How- 
ever, these discrepancies do not invalid the usefulness of 
the ten-spin model for appreciating the effects on the spin 
dynamics of the many-body quantum character of the 
spin interactions, like the eigen-selectivity, as was shown 
in Section II A. 

The measurements in PAA^e and PAA were made at 
temperatures T = 115°C and T = llCC, respectively. 
In the remaining cases, the experiments were carried out 
at room temperature, namely T = 27°C. The MREV8 
sequence shown in Figure 1 (b) was set up to mitigate 
the effects of the finite width of the pulses. The total 
time of each sequence is Tc = 2 (2r2 + 2ti + 4t^)^^, then 
T\ = Tc/l2 — tw y T2 = Tcl^ — ty^. For a pulse of 7r/2 with 
a width of = 6.56/is and a minimum setting time of 
Ti — 1.6^s, we have Tc — 12 (ti +t^) — 97.92fis and T2 — 
9.76/xs. By comparing the time scale of the FID's shown 
in Figure 5, it can be anticipated that the evolution under 
the dipolar Hamiltonian during the pulse duration will 
not have relevance in the experimental results. 

In Figure 6 we present experimental results corre- 
sponding to the refocalization experiments discussed in 
Section II A in 5CB under the sequence shown in Figure 1 
(c). From the spectral evolution of the 'time-reversed' 
FID's of Figure 6 (a), given in Figure 6 (b) or in more 
detail in Figure 6 (c) and Figure 6 (d) . it can be seen how 
the spin dynamics is affected by the eigen-selection pro- 
cess during the time reversion (in the last figure the spec- 
tra as a function of t are normalized to facilitate the com- 
parison). This becomes more evident when comparing 
the amplitude variation of two groups of spectral lines, 
those around the frequencies 5.55KIIz and 10.60KHz, as 
shown in Figure 8 (a). There it can be seen that the peak 
at 10.60KHz presents a higher decay rate for all t. It can 
be appreciated in Figure 6 (c) how the two peaks match 
their amplitudes at r « 361.4/is, reverting subsequently 
their initial amplitude relation. It is worth to note that 
in the spectra the frequencies lines close to f=OHz are 
affected by instrumental artifacts like noise in the base- 
line of the signals. Therefore, the variations in such line 
should not be giving any physical meaning. 

The results of the same experiment performed in ne- 
matics FAA^e and PAA, and solid adamantane, are 
shown in Figure 7. In PAA^^g, in agreement with the re- 
sults of the numerical calculation shown in Section II A, 
it is observed a narrowing of the spectrum, in consistency 
with a more pronounced decay of the higher frequencies 
lines. By other hand, in the PAA spectrum low and 
high frequency groups of lines are resolved. Figure 8 (b) 
shows the detail of the normalized decay of two lines at 
l.GOKHz and 7.35KIIz. It should be noticed that the de- 
cay rate depends on the line position in the frequency 
spectrum, and the 'bell' form associated with the homo- 
geneous broadening of every line (the lincshapc produced 
by AQD) is equally affected by decoherence or non-ideal 



effects like static field inhomogeinity. Due to this, for gen- 
eral cases the time evolution of the different parts of the 
spectrum will be masked by the superposition of broad- 
ened lines, preventing the experimental resolution of the 
decay rates of different parts of the spectrum. Therefore, 
it can be expected to obtain a suitable detection of the 
eigen-selectivityin samples with spectrum having groups 
of lines with certain degree of resolution, like nematics 
5CB and PAA. 

The experiment sketcehd in Figure 1 (b) can be mod- 
ified by changing the MREV8-train by a single block, 
where the time parameter ti is varied independently, and 
the value of T2 is defined as to obtain the maximum sig- 
nal, being T2 > n . For these values of T2 as a function of 
Ti the slope of a stright line is fitted, whose ideal value 
is two. By using this dependence, the reversion exper- 
iment can be done by using a single block, in such a 
way to minimize non-idealities present in the pulse set- 
ting of the experiments. In this experiment, the rever- 
sion time r is a linear function of ti. The advantage 
of this method is that it permits a continuous variation 
of Ti yielding a smooth profile for the signal amplitud 
in the reversion experiments. On the other hand, long- 
time settings in such experiment can introduce contri- 
butions from the non-secular dipolar Hamiltonian on the 
dynamics (see supplementary matcrial^^, section "Misad- 
justments in the dynamics under reversion", for details 
about the influence of the non-secular dipolar part on the 
MREV8 reversion dynamics). 

The results of the reversion experiment using the con- 
tinuos MREV8 method in 5CB are shown in Figure 9. In 
Figure 10 is shown the normalized variation with the re- 
version time of the amplitude of two peaks at distinct fre- 
quencies, corresponding to the measurement of Figure 9. 
There it can be appreciated the frequency selective de- 
cay. Due to the finite pulse widths, the initial measured 
amplitudes correspond to a time r significantly shifted 
from zero. For these reversion times the frequency se- 
lective decay is noticeably, accordingly the first spectra 
in Figure 9 already presents a larger decay of the higher 
frequency. This explains the difference with the spec- 
tra shown in Figure 6. Except for this last detail, the 
obtained results show the same features than the rever- 
sion experiments using blocks of the MREV8 sequence. 
Similar results were obtained for the other compounds 
studied in this work. 

By means of the continuous MREV8 experiment it is 
possible to obtain a detail of the maximum amplitude of 
the reversed FID's, as a function of the reversion time r. 
In Figure 11, the results obtained by averaging data near 
the maximum of each reversed FID are shown, where the 
signal to noise ratio is larger. Besides, the plots were nor- 
malized with respect to the first FID obtained in each ex- 
periment. We added three gaussian profiles with different 
values of standard deviation {a = 350/Lts, 580/xs, 800/zs), 
used as a guide to the eye. Since the first FID signal 
already present a significant attenuation, the gaussian 
curves have a value greater than one for r = 0. It can be 
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seen that the decay time is longer than the characteristic 
decay time of the FID for every compound (see Figure 5). 
This is consistent with the theoretical approach proposed 
in QD-I, where it is considered that the spin dynamics 
during the reversion period is governed by a different 
mechanism than the one which control the spin dynam- 
ics during the FID evolution. This quantum process is 
characterized by a longer time scale, since it is associated 
with higher order terms of a perturbative treatment. Be- 
sides, we can see that for times shorter than 780/is the 
PAArfe (with 8 spins per molecule) presents a similar de- 
cay behavior that the PAA (with 14 spins per molecule) 
and the 5CB (with 19 spins per molecule) presents a sim- 
ilar decay behaviour that the adamantane (a solid array 
of spins) for all r time. These results reflect that deco- 
herence is not associated with the nature of the spins as 
a closed system. 

It is also observed in Figure 11 that the decay in PAA^e 
is more similar to that of PAA than to the one of 5CB. 
The fact of observing similar responses to decoherence 
for these two samples is consistent with decoherence be- 
ing controlled by the coupling of the spin system to ex- 
ternal degrees of freedom. For values of t greater than 
780/Lts, decoherence in PAA turns stronger, which could 
be indicating the occurrence of an additional decoher- 
ent mechanism associated with the protons of the methyl 
groups. However, because of the low signal to noise ratio 
of PAAdG and other sources of error, which will be dis- 
cussed in Section II C, more experiments are necessary to 
confirm this assertion. 

It is worth to note that, using a single MREV8 pulse 
sequence with a continuous variation of time settings, 
for long r values the decay can be affected by the 
influence of the non-secular dipolar Hamiltonian part 
in the dynamics under reversion (see supplementary 
materiaP'', section "Misadjustments in the dynamics 
under reversion", for details about the influence of such 
non-secular dipolar terms on the MREV8 reversion 
dynamics). However, such experiment is able to show 
the intermediate time scale of decoherence as it is the 
aim in this work. For a more exactly measurement of 
the decoherence time other techniques can be used, like 
the 'magic sandwich'^'^'^^'^^""^'', when they are adequate 
due to constraints in the experimental setting like the 
duration of the rf pulses. 

Summarizing, from the results presented in this sec- 
tion, we have the following conclusions: 

1. the time scale of the attenuation of the time re- 
versed FID signal in MREV8 experiments occurs 
in a longer time scale than the one of the FID evo- 
lution, 

2. the high-frequency spectral components decay 
faster, showing an eigen-selection process, 

in total agreement with the theory presented in QD-I. 
Besides, such intermediate time scale for decoherence. 



between the FID evolution and the thermalization 
process time scales, is consistent with the hypothesis 
about the existence of dynamics with different time 
scales used in the theoretical approach of this work and 
QD-I. 



C. Non-ideal behaviors and approximations in the 
experiments 

The MREV8 sequence shown in Figure 1 (b) has been 
optimized by adjusting the time parameters ti and T2, to 
mitigate the effects of the finite width tyj of the pulses, 
as was commented in Section II B. The shortest interval 
Ti between pulses was fixed in 1.6/zs. It is worth to note 
that in a time of 6.56/iS, corresponding with a 7r/2 pulse 
in our experiment, the free evolution dynamics produced 
by the dipolar interaction is not significant as can be seen 
from the FID's shown in Figure 5. 

The non-ideality arising in the control of the reversion 
parameter k (see Section II A for details about this re- 
version parameter) would come from deviations from the 
exact value k=2 due to misadjustments of the 7r/2 pulses. 
In a similar fashion, a correct setting of the tt/2 pulses 
with K=2 could be affected by an incorrect adjustment 
of the time T2 (i.e. with T2 ^ «ti), causing an imperfect 
reversion of the spin dynamics. These errors could intro- 
duce an additional dynamics of the 'essentially closed' 
system^^, generated by the molecular dipolar Hamilto- 
nian, producing perturbations of the signals and their 
spectra as a function of the reversion parameter r. 

With the aim of estimating the effects of experimental 
misadjustments, in the following we will obtain several 
analytical expressions for different kind of errors in the 
reversion sequence shown in Figure 1 (c) using MREV8 
blocks of pulses. We will consider the spin system as a 
closed system to study the possibility of obtaining some 
eigen-selection effect in such case, which could overlap 
with the effects of the coupling with the environment. We 
relegate to the supplementary material^^, section "Mis- 
adjustments in the dynamics under reversion" , the com- 
plete analytical demonstration, as well as the numerical 
calculation of the signal expressions which will be used 
in this section. 

We introduce the error e in the setting of the time un- 
der the reversion dynamics by writting t2 = {k + e) n. 
It is worth to note that such error takes account of mis- 
adjustments in the time T2 as well as deviations from 
a perfect setting of the tt/2 pulses (see supplementary 
material^^, section "Misadjustments in the dynamics un- 
der reversion" , for details about the equivalence between 
T2 time misadjustments and 7r/2 pulses misadjustments). 
Therefore, we have ^2 = (^ -I- e) ti and the total reversion 
time r = (K-|-l-|-e) ti. Using such erroneous setting of the 
t2 time, we obtain for the expectation value of ly corre- 
sponding to the closed spin system, under 'on-resonance' 
condition, the following result 
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where we used the approximation c/k <C 1 which imphes 
t\ '^t/{k + 1). The Fourier transform of (71) is 
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The expression in Eq. (72), for the Fourier transform of 

the FID signal for a closed spin system, does not present 
a decay with some eigen-selection effect, but it shows an 
oscillatory behaviour of each spectral line proportional to 
the product of the eigenvalue difference — Ci with the 
error factor e/n. 

Other interesting case; to consider, is when the FID signal 
is considered as produced by the different signals com- 
ing from a distributed molecular orientation of the main 
molecular axis of the LC system'*-'^. In such case, with 
the same misadjustment as in Eq. (71), we have for the 
FID 
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where p\^'' (Si) is the OMDF (see Section II A). The ex- 
pression (73) is the classical representation of the AQD 
seen in Section II A. Since the dynamics produced by 
the error factor c/k under the reversion time r is very 
slow in comparison with the one produced under the free- 
evolution time t, we can write 
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The Fourier transform of (74) is 
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We can see from Eq. (75) that there is a decay factor 
with eigen-selection effect for each spectral line, but it 
is masked by an oscillatory behaviour of each frequency 
line, instead of the monotonous decay that is experimen- 
tally observed (see for instance Figures 6 and 8) . Besides, 
such predicted decay is reversible and it would, in princi- 
ple, be possible to increase the decay time by correcting 
the misadjustment represented by the factor e/n, since 
the dynamics is very susceptible to changes of e. Such 
changes of e modify the decay rates for each spectral lines, 
but the oscillatory behaviour of such lines can not be sup- 
pressed in this way. Indeed, this effect is not observed in 
the experimental measurements, where the correction of 
the misadjustments produces the extinction of the oscil- 
latory behaviour under the reversion dynamics. 
If we see the signal under the reversion time r obtained 
for the time t = Q, we have 
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The expression (76) corresponds to an FID signal but 
with a slower time behaviour in comparison with the free- 
evolution FID. Such behaviour under reversion depends 
on the time r instead of the time t. The dynamics 
represented by Eq. (76) is reversible and the oscillations 
in the signal are unavoidable. Clearly, this is not a fea- 
ture of the experiments, as was mentioned. 
If the setting of the 7r/2 pulses is very erroneous or if the 
dynamics under reversion is very influenced by the ef- 
fect of the non-secular dipolar Hamiltonian (for instance, 
when the time ti is not small enough to neglect the 
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dynamics produced by the non-secular dipolar Hamilto- 
nian), we can obtain a dynamics for the signal under 
reversion as 
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where ^5*^' is the resulting Hamiltonian under re- 
version, which is different of H^gi"^ (see supplementary 
material'^^, section "Misadjustments in the dynamics un- 
der reversion", for details of the definition of H^g^^^). In 
Eq. (77) the eigenbase { | ai )} of the Hamiltonian H^si^^ 

is defined, where ^si^'jcn) = QJijcKi)- The Fourier 
transfor of (77) is 
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where we can see that the decay corresponding to the 
spectral lines does not present any eigen-selection effect. 
These effects of misadjustments in the reversion time, 
with the spin system considered as a closed system, 
were numerically simulated and the results are shown in 
the supplementary material^^, section "Misadjustments 
in the dynamics under reversion", where we obtain 
numerical conclusions which coincide with the analytical 
analysis. 

Finally, to contrast with the commented analytical re- 
sults about misadjustments in the experimental setting 
in the closed spin system, we derived the analytical signal 
produced by the same reversion experiment affected by 
an erroneous setting of the reversion time, for the open 
quantum system conformed by the protons coupled to 
the environment. Using the time setting t2 = {k + e) ti 
in the evolution operator (43), we obtain the decoherence 
function 



(79) 

instead of the Eq. (48), where we have used the same 
considerations involved in Eq. (48). In Eq. (79), we de- 
fined 



yAs{t,Ty = t 



t/ K, 



yL{t,Ty = T^iy{f + rf-[l + {K + l + e)-^] ■ 
-[{K+1 + e)-i + (« + 1 + e)-2] i T^j, 



ysL{t,'^ 



(80a) 



(80b) 



(80c) 



Under the approximation in the dynamics of 
/// riz) rt {■}t=0' used to obtain the well approxi- 
mated decoherence function shown in Eq. (53), and with 
c/ac <C 1, we can write (79) as 

G,„.,((,r)'^ / <iAS.e-«.-«>^*(-*') 

(81) 

The decoherence function (81), affected by a misadjust- 
ment of the reversion time will affect the expectation 
value of ly in the way shown in Eq. (28). We can 
see from Eq. (81) that the decays produced by the 
integration of the complex exponential function under 
the total reversion time r can be compensated excepting 
for the dynamics under the eigenvalue Cf' . Such decays 
present the eigen-selection effect and the oscillations 
in the signals under the reversion time r, due to the 
misadjusment of ^2, can be cancelled. It can also be 
seen that the decays associated to the eigenvalues A5i 
and Cfl^ can be reverted but the decay due to will 
never be reverted. Such non-reverted decay constitutes 
an enveloped which constrains the signal as a function 
of the reversion time t giving it a bell-like shape, as we 
will shown in the following text about the experimental 
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measurements. By making e = in Eq. (81) we recover 
Eq. (53), as expected. 

In Figure 12 (a) we show the effect mentioned above for 
the reversion experiment in 5CB, where in Figure 12 
(b) it can be seen the variations of the amplitudes of 
the pseudo-FID, obtained extracting the mean values of 
the amplitudes of the FID around t = 0, introduced by 
the error of k. By comparing Figure 12 and Figure 6, 
it can be noticed that in Figure 6 (a) the amplitudes 
decay monotonously, indicating that the error k has 
been satisfactorily mitigated. This dynamics induced 
by experimental mismatches cannot produce a definitive 
signal decay, because of the small number of spin degrees 
of freedom. The latter is true for any error in the pulse 
configuration. In practice, this effect can be attenuated 
by varying the time between the two WHH-4 blocks 
comprising the MREV8 sequence until the oscillations of 
the amplitudes as a function of r dissapear (see Figure 1 
(b)). The results shown in Figures 6 and 7 were obtained 
with this procedure. 

The observed signal attenuation might be associated 
with several causes besides the essentially adiabatic 
decoherence^^, like fluctuations of the spin-spin interac- 
tions due to thermal molecular motions^^ or even experi- 
mental non-idealities like inhomogeinity of the static and 
the rf magnetic fields. In the work of reference^"^, the ex- 
perimental sources of error in the 'magic-sandwich' (MS) 
reversion experiment were carefully checked by study- 
ing the effect of the sequence on the spin system in the 
isotropic phase of liquid crystal 5CB. Combination of the 
MS sequence and a tt pulse to reverse the static field in- 
homogeinity yielded the same response than the usual 
Halm-echo two pulse sequence (T2 = 70ms). This also 
showed that the time scale of the decay produced by ther- 
mal fluctuations is much greater than that of the MS 
experiment in the ncmatic phase. 

In the present work, we studied analytically and ex- 
perimentally the influence of the field inhomogeneity in 
the FID dynamics and the reversion dynamics perform- 
ing experiments on isotropic 5CB (see supplementary 
material^^, section "Effects of the field inhomogeneity", 
for details about such experiment). The experimental 
results allowed us to conclude that the dynamics pro- 
duced by the field inhomogeneity has a time scale longer 
than that of decoherence and it does not present a be- 
haviour with eigen-selectivity. Besides, it is mentioned 
that MREV8 experiments combined with tt pulses were 
carried out by the authors in 5CB in the isotropic phase 
where it was also obtained an exponential decay with a 
T2 very similar to the obtained in 5CB in reference 
On the other hand, the homogeinity of the rf field was 
optimized by using a low coil filling factor and the signals 
obtained from samples with different sizes do not present 
any different behaviour between them. 

These observations give support to the statement that 

while 'bulk' effects like inhomogeneity of the rf pulse and 
thermal fluctuations of dipole-dipole interaction could 



eventually perturb the free spin dynamics, their influ- 
ence is irrelevant in the intermediate time scale of the re- 
version experiments. Therefore, the observed decay can- 
not be associated with experimental non-idealities, but it 
should be assigned to the irreversible spin dynamics in- 
duced by quantum decoherence, whose fingerprint is the 
eigen-selectivity. 

Summaryzing, even when the observed signals are af- 
fected by magnetic field inhomogencities and intcmolec- 
ular dipole interactions, the contribution of the quantum 
molecular mechanical dynamics is more important along 
the different time scales. 



III. DISCUSSION AND CONCLUSIONS 

In this work, we studied the spin dynamics which 
characterizes the irreversible decoherence of a finite 
quantum interacting spin system coupled with an infi- 
nite quantum environment. We experimentally detected 
in nematic liquid crystals the salient characteristics of 
the spin dynamics predicted by a full theory which 
considers the spins as an open quantum system, namely 
eigen-selectivity, spectral compression and irreversible 
decoherence under refocusing of the dipolar spin inter- 
actions. 

The experiments were interpreted in the context of 
the theory presented in the work of reference"'^"'^, under 
the main assumption that irreversible decoherence 
occurs well before that thermal fiuctuations play any 
significant role, and long after the dephasing by quantum 
interference and reversible adiabatic decoherence. This 
hypothesis is supported both by our experiments and 
the work of reference-'^^. where it was shown that irre- 
versible decoherence, which cannot be associated with 
thermalization (neither adiabatic nor nonadiabatic), 
occurs in an intermediate time scale. Besides, the 
accurate description achieved in reference^^ of the time 
domain FID signal by including quantum interference 
and reversible adiabatic decoherence, allowed us to 
show that the time scale characterizing the Liouvillian 
dynamics and such reversible decoherence (which can 
be interpreted semiclassically) is much shorter than the 
time scale of the irreversible quantum decoherence. 

The analytical- numerical treatment of Section II A 
allowed us to compare the theoretical expressions derived 
from the proposed theory, with the experiments. By 
introducing decoherence functions with gaussian decay 
profiles, in the free-evolution and reversion dynamics, the 
effect of the environment was included. The calculations 
were performed on a 10-spin model for 5CB^^ and on 
the 8-spin model of PAA^g presented in QD-I-'^^. The 
results of the calculations are shown in Section II A and 
the measurements are shown in Section II B. By using 
reversion experiments we confirmed the occurrence of the 
intermediate time scale and the characteristic behaviour 
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of the irreversible decohcrcncc, and comparison of the 
analytical FID signals and their spectra under reversion 
dynamics with the experimental ones confirmed the 
validity of the theoretical approach presented in QD-I^^ 
and in this work, beyond the experimental non-idealities. 

We present a detailed analysis of the experimental 
causes that could affect our measurements, in order 
to identify their characteristic time scales. In the 
supplementary materiaP^, section "Effects of the field 
inhomogeneity" , we analyse the effect of the inhomo- 
gcneitics of the static magnetic field. We conclude that 
the dynamics produced by this effect has a longer time 
scale than decoherence, and most importantly, it does 
not present eigen-selectivity. With the aim of checking 
for a possible frequency dependent behavior induced 
by pulse misadjustmcnt that might be confused with 
genuine eigen-selectivity, in Section II C we present a 
theoretical analysis of the effects of these misadjustments 
on the FID signals, considering the spins as a closed 
system. The simulations on PAA^e presented in the 
supplementary material, section "Misadjustments in the 
dynamics under reversion" show that misadjustmcnt 
cannot be the source of the observed eigen-selectivity of 
a closed spin system. Besides, this kind of non-idealities 
were theoretically analyzed for decoherence produced 
by the full-quantum dynamics in Section II C, as well. 
The behaviour of the signals extracted by such analysis 
agrees with the measurement observed. Therefore, this 
detailed analysis leads us to conclude that the observed 
eigen-selectivity is an evidence of the open quantum 
system dynamics. 

Eigen-selectivity, which introduces a distinction in 
the response of the diagonal and off-diagonal elements 
of the density matrix, provides an efficient irreversible 
mechanism for coherence decay, while preserves the 
'population' terms which can only change in a much 
longer time scale. This behavior explains the observed 
buildup of the quasi-equilibrium in liquid crystals. 
The results of this work, together with the conclusions 
obtained in references^^'^'\ can contribute to elucidate 
the quantum mechanisms underlying decoherence of 
open quantum systems of interacting spins, for instance 



the role played by quantum correlations between the 
observed system and the environment in the damping 
of the spin cohercnces^'''^^. Certainly, these works"^""^'^^ 
show that in liquid crystals it is essential to assume the 
quantum character of the spin-environment coupling to 
explain the observed irreversible decohercince. These 
statements might also apply to interacting spins in 
ordinary solids, where irreversible adiabatic decoherence 
could provide an explanation for the quasi-equilibrium 
states characterized by spin temperatures^^. 

Summarizing, the occurrence of eigen-selectivity in the 
spin dynamics under reversion was verified in several ne- 
matic liquid crystals, through the direct experimental 
observation of inhomogeneous decay and spectral com- 
pression of the NMR spectrum under refocusing of the 
dipolar spin interactions, over an intermediate time scale. 
We conclude that the eigen-selection effect is the finger- 
print of decoherence associated with a quantum open 
spin system in liquid crystals. Besides, the dynamics of 
such interacting spins, with few degrees of freedom, can 
be described in terms of quasi-equilibrium states of each 
molecule after the irreversible decoherence damps out the 
coherent part of the spin state. Therefore, the observed 
system reaches these states through a genuine quantum 
process involving spins and environment, being this a 
process of different nature than the fluctuations which 
govern thermalization and relaxation. These findings al- 
low to understand the buildup of the quasi-equilibrium 
states as being a consequence of the correlated dynamics 
of the observed system and the quantum environment. 
Accordingly, the quasi-equilibrium representation in liq- 
uid crystal needs not being perceived heuristically, in- 
stead, it should be considered as a definite stage of the 
spin system, during its evolution towards equilibrium. 
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FIG. 1: Experimental pulse sequences, a: Single reversion, 
b: MREV8. c: Compound reversion sequence, where each 
reversion block RS{i) constitutes a single reversion sequence 
like the MREV8. 




FIG. 2: Calculation of the FID signals under spin rever- 
sion dynamics, using MREV8 pulse blocks, for the nematic 
5CB molecule in the resonance condition. Decoherence is 
represented with gaussian profiles in t and r. Parameters: 
S^^ = 0.5403, at = 0.07 (standard deviation in t), cr^ = 0.04 
(standard deviation in t). a: FID signals, depending on time 
t, as a function of the reversion time r. b: Amplitude spectra 
for the results in (a), c: Frontal detail of the spectra shown in 
(b). d: Spectra of (c) normalized to compare the evolution of 
their frequency components (for clarity only the spectra for 
the shorter r values are shown). It is observed in (d) the dif- 
ferent decay rates of a set of spectral lines close to 5.65KHz 
and 8.50KHz as well as the spectral compression which are 
evidences of the eigen-selectivity in the decoherence process. 
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FIG. 3: Calculated evolution of the normalized amplitude as 
a function of the reversion time r for the spectral lines at 
5.65KHz and 8.50KHz of the spectra shown in Figure 2 (b) 
and (c) for the nematic 5CB. 




FIG. 4: Calculation of the FID signals under spin rever- 
sion dynamics, using MREV8 pulse blocks, for the nematic 
PAA^e molecule in the resonance condition. Decoherence is 
represented with gaussian profiles in t and r. Parameters: 
Szz ~ 0.53, at = 0.06 (standard deviation in t), (Jt ~ 0.015 
(standard deviation in r). a: Amplitude spectra for the cal- 
culated FID's. b: Spectra of (a) normalized to compare the 
evolution of their frequency components (for clarity only the 
spectra for the shorter r values are shown). It is observed in 
(b) the spectral compression which is evidence of the eigen- 
selectivity in the decoherence process. 
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FIG. 5: Molecular sketch (inner detail in the figures), FID 
evolution (index 1 in the figures) and the FID spectrum (in- 
dex 2 in the figures) measured for nematic 5CB (al,a2) at 
room temperature (namely T = 27°C), PAAde (bl,b2) at T 
= 115°C, PAA (cl,c2) at T = 110°C, and solid adamantane 
(dl,d2). 




FIG. 6: Experimental results of the FID signals under spin 
reversion dynamics, using MREV8 pulse blocks, for the ne- 
matic 5CB at room temperature (namely T — 27°C) in the 
resonance condition, a: Measured signals for the FID experi- 
ments, depending of the time t, as a function of the reversion 
time T. h: Amplitude spectra for the measurements shown in 
(a), c: Frontal detail of the spectra shown in (b). d: Spectra 
of (c) normalized to compare the evolution of their frequency 
components (for clarity only the spectra for the shorter r val- 
ues are shown). In (c) and (d) is detailed the time r for each 
spectrum. 
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FIG. 7: Experimental results of the FID signals under spin 
reversion dynamics, using MREV8 pulse blocks, in the reso- 
nance condition. The figures with index 1 show the amplitude 
spectra for the FID measurements and the ones with index 2 
show the measured FID spectra, which are normalized to com- 
pare the evolution of their frequency components (for clarity 
only the spectra for the shorter r values are shown). al,a2: 
PAAds at T=115°C. bl,b2: PAA at T=110°C. cl,c2: Solid 
adamantane. 
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FIG. 8: Measured evolution of the normalized amplitude as a 
function of the reversion time r for different spectral lines of 
the FID spectra, a: Spectral lines at 5.55KHz and 10.60KHz 
for the spectra shown in Figures 6 (b) and (c) for nematic 5CB 
at room temperature (namely T — 27°C). b: Spectral lines at 
l.GOKHz and 7.35KHz for the spectra shown in Figure 7 (bl) 
for the PAA at T = 110°C. 




FIG. 9: Experimental results of the FID signals under spin re- 
version dynamics, using a single MREV8 pulse sequence with 
a continuous variation of its setting times, for nematic 5CB 
at room temperature (namely T — 27°C) in the resonance 
condition, a: Measured signals for the FID experiments, de- 
pending of the time t, as a function of the reversion time r. 
b: Amplitude spectra for the measurements shown in (a), c: 
Frontal detail of the spectra shown in (b). d: Spectra of (c) 
normalized to compare the evolution of their frequency com- 
ponents (for clarity only the spectra for the shorter r values 
are shown). 




FIG. 10: Measured evolution of the normalized amplitude as 
a function of the reversion time r for the spectral lines at 
5.60KHz and 8.35KHz for the spectra shown in Figures 9 (b) 
and (c) for the nematic 5CB at room temperature (namely T 
— 27°C). It can be seen the eigen-selectivity where the higher 
the frequency of the line is, the faster its decay will be. The 
continuous variation of the setting time allows to see clearly 
the bell-like form for the decay of the spectral lines. 




FIG. 11: Experimental results of the normalized evolution of 
the maximum of the FID signals under spin reversion dynam- 
ics, using a single MREV8 pulse sequence with a continuous 
variation of its setting times, for adamantane and 5CB at 
room temperature (namely T = 27°C), PAAds at T = 115°C 
and PA A at T = 110°C. Three gaussian decays with different 
standard deviation, a, are plotted for comparison with the 
measurements. 




2500 
3000 Tins] 



b 1'- 

0.9 \ 









1 : 1 1 ; 
1 ; 1 1 ; 






1 : : 1 ; 
1 : : : : 

1 ; ^ 






1 : i i ; 

--t -J i -; 






1 : : 
1 : 1 : 






^ i 1 I ; 






— \ - --- - 

\ . 

\ : : 






V : 1 : 






"k \ ; i 






: i ; 






: V :...jKrT.j L 




1 1 1 1 1 





1 1 I \ \ I I 1 I 

500 1000 1500 2000 2500 3000 3500 



FIG. 12: Experimental results with a misadjustment in the 
reversion parameter k of the FID signals under spin rever- 
sion dynamics, using MREV8 pulse blocks, in the resonance 
condition for nematic 5CB at room temperature (namely T 
— 27°C). a: Measured signals for the FID experiments, de- 
pending of the time t, as a function of the reversion time r. b: 
Pseudo-FID obtained for the evolutions in r of the FID-signal 
values in (a) at t — 18/is. Such pseudo-FID presents a signal 
form similar to the FID signal, but with lower characteristic 
frequencies. 
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The material below is intended as supplementary 
material to the article Quantum irreversible de- 
coherence behaviour in open quantum systems 
with few degrees of freedom. Application to 
NMR reversion experiments in nematic liquid 
crystals. Such manuscript will be referred to as QD-II 
in the following text. 

In Section I we analytically and experimentally will 
analyse the time scale of the effects produced by the in- 
homogeneity of the strong magnetic field on the signals 
obtained from spin obscrvables in liquid crystals (LC). 
Section II is dedicated to estimate the effects of mis- 
adjustments of the characteristic times of the reversion 
pulse sequence, as well as the influence of the non-secular 
parts of the dipolar Hamiltonian in such experiments. 

I. EFFECTS OF THE FIELD INHOMOGENEITY 

This section is concerned with the analysis of the field 
inhomogencity effects on the experimental measurements 
presented in QD-II. The field inhomogencity causes the 
decay of the NMR signals due to a distribution of the 
z component of the static magnetic field. Therefore, it 
is important to know the characteristic behaviour and 
time scale of that decay in order to compare with the 
decay produced by decoherence. 

To observe the effects of the field inhomogeneity in the 
nematic phase, we first study the FID signals obtained 
by applying the Zceman reversion sequence shown in Fig- 
ure 1 (a). In a nematic LC, the static magnetic field docs 
not vary significantly inside of each molecule but there 
is a distribution of the field strength along the whole 
sample. Such distribution produces a dispersion of the 
Larmor frequency in the sample. Therefore, the evolu- 
tion operator due to the field inhomogeneity for a single 
molecule in the 'on resonance' condition will be 

U^(t) = e-^^^-*, (f ) 

where A = co — ojq. The frequencies w = 7B and 
oJo = jBq, with the proton gyromagnetic ratio 7, are 

produced by the strength B of the magnetic field inside 
the molecule and the mean value Bq of the strength of 



the magnetic field of the sample, thus wq is the frequency 
of the rotating system and it defines the resonance con- 
dition. 

The complete dynamics in this condition for the inter- 
mediate time scale (see QD-II for details about the ex- 
istence of different time scales in the dynamics) will be 
produced by the secular full-quantum dipolar Hamilto- 
nian and the environment or lattice Hamiltonian, stud- 
ied in QD-II, together with the Zeeman-likc Hamiltonian 
in Eq. (1). Since the Hamiltonian in the evolution op- 
erator (1) commutes with the other ones, we can study 
the dynamics produced by such evolution operator inde- 
pendently of the other evolution operators and readily to 
include the effects of the inhomogeneity on the complete 
dynamics. 

Keeping the last statement in mind, we can extract the 
evolution operator due to the field inhomogeneity of the 
dynamics under the Zeeman reversion experiment shown 
in Figure 1 (a) for a single molecule and then to conclude 
how it affects the complete dynamics. Such operator is 
written as 

U|r(r) = R-]." e-'^^-^/2 e-'^^-^/2 = 1, (2) 

where is defined the operator = e* of the pulse in 

an angle tt and direction x and we used Rj^IzR-x = ~Iz- 
Therefore, the expression in Eq. (2) tells us that the 
field inhomogeneity effect in the whole sample will be 
removed and the complete dynamics will produce the 
corresponding signals for the observables, for instance 
the FID, in the same way as if the inhomogeneity would 
not exist for an ideal experiment. 

We can see the results of such experiments in Figure 2 
for the nematic LC's 5CB and FAA^e- As is observed 
there, and comparing with the measured FID's shown in 
Figure 5 in QD-II, there is no change in the time scale 
of the FID signals under Zeeman reversion. Therefore, 
there is not an influence of the field inhomogencity 
whose effects can be observed in the FID time scale. It 
is worth to note that the FID's observed in Figure 2, 
under a Zeeman reversion dynamics, are different than 
the free-evolution FID's. Such differences arc produced 
by experimental issue, like the finite width of the tt 
pulses, but they do not change the time scale of the 
output signals nor the conclusions extracted about that. 
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Then, wc studied the effect of the field inhomogcneity 
over the dynamics of the pseudo-FID's obtained under 
the reversion pulse sequence shown in Figure 1 (c) using 
as RS{i) blocks an MREV8-like sequence as shown in 
Figure 1 (b). The aim of such reversion experiment, in 
our work, is to remove the dynamics produced by the 
dipolar Hamiltonian. Nevertheless, the environmental 
decoherence, which is studied in QD-II, and the non- 
ideal effects or misadjustments in the reversion dynam- 
ics, which are studied in Section II in this supplementary 
material, avoid the complete removal of the dipolar dy- 
namics and a final decay of the signals of the observables 
is measured. The field inhomogcneity will introduce an 
extra dynamics produced by a Zceman-like Hamiltonian 
(1), and that dynamics will be observed as an extra de- 
cay in the measured signals (as we will see), which will 
be mixed with the other mentioned decays. 

To extract the time scale of the decay due to the field 
inhomogcneity under the mentioned reversion sequence, 
we performed the experiment described in Figure 1 (c) 
in a 5CB sample in isotropic phase. In such phase, the 
dynamics is governed by a Zeeman-like Hamiltonian of 
the inhomogcneity. which is the same Hamiltonian as in 
the nematic phase, but in the isotropic phase the dipo- 
lar Hamiltonian is averaged out. In the longest time 
scale, the thermalization, which is studied by relaxation 
theories, will produce a final decay of the signals, but 
such time scale is much longer than that of decoherence. 
Therefore, the fast decays observed in the pseudo-FID's 
under reversion dynamics in the isotropic phase can be 
associated meanly to the field inhomogcneity. 

In the same way that in the former experiment, where 
it is obtained the evolution operator (2), we first calculate 
the evolution operator due to the field inhomogcneity for 
the sequence shown in Figure 1 (b) for a single molecule. 
Such evolution operator is written as 



signals due to the field inhomogcneity. 

Next, using the evolution operators for a free-evolution 
dynamics, Eq. (1), and for a reversion dynamics, Eq. (4), 
we obtain the signals of the whole sample of an isotropic 
LC affected by a field inhomogcneity as 

X g±(CSi,Ci4)»A.(t+r/3) 

(5) 

where the index i labels the i-th molecule (the sum runs 
over the molecules of the whole sample), {| G^i )} is a 
base in the Hilbert space of the spins of the i-th molecule 
which are eigenstates of the spin part of the interaction 
dipolar Hamiltonian. Since the dipolar and the Zeeman 
Hamiltonians commute, such base can be selected as a 
common eigenbase of the Zeeman Hamiltonian as well. In 
Eq. (5), the symbol ±{Qs^,C,iSi) represents the possible 
values or —1 depending of the difference in the Zee- 
man eigenvalues of the states | Ci^i ) and | Ci^i )' ^^^^ 
ference can only have the values ±1 due that the non-null 

(s ■ ) 

matrix elements of the operator ly ^' connect eigenstates 
with such difference of Zeeman eigenvalues. Besides, we 
defined Afsi = trg^ {l^^^^} and (St = -j^^, with ks the 
Boltzmann constant and T the absolute temperature. 
The matrix elements ( Qs^ |ly '''I Ci-^i ) are the same for 
different molecules. We suppose a dense and continu- 
ous distribution of the Larmor frequency shifts A for the 
molecules (which is realistic for a typical field inhomo- 
geneity). Due to this, we can replace the sum over the 
molecules in the sample as an integration over the values 
of A. Therefore, Eq. (5) can be written as 



TJ^ c(ti) = p-*^!^"^! p^'^Iy^i „2AI^2ri zAIyTi 

^ g-iAI^2Ti g-iAIyTi g-iAI^2ri g-iAIyn g-iAI^n 

It can be observed in Eq. (3) that there are exponen- 
tial operators with the same exponent and with opposite 
sign, thus it will produce a cancellation effect of their 
dynamics. For small enough values of the time ri, the 
operator (3) can be approximated as 

Uis(r) ^ e-'^^^-/^ (4) 

where r = 12ri for a single MREV8 block. A reversion 
experiment with n blocks of the MREV8 pulse sequence 
will have for the field inhomogcneity a similar evolution 
operator as in Eq. (4), in the approximation of small 
times, with r = 12nTi. 

Comparing Eq. (4) with Eq. (1), we can see that the 
dynamics under reversion using MREV8 blocks is a 
third slower than the free-evolution dynamics, thus the 
reversion experiment mitigates the decay effect over the 



X Gz(t + r/3). 

(6) 

In Eq. (6) it is defined 

Gzit)= J dAe^^Vz(A), (7) 

where pz{A.) is the distribution function of A and 
pz{A)dA — n{A)/N with n(A) as the number of 
molecules with the value of A and N as the number 
of molecules of the whole sample. Given that pz(A) is 
supposed a symmetric bell-form distribution with a null 
mean value (which is realistic as well), the integration 
(7), whose limits of integration can be extended to ±oo, 
produces the same result if the complex exponential has 
positive or negative sign, thus it is independent of the val- 
ues of ±{(iSi,([s[) and we only keep the positive sign. 
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To extract an expression for Eq. (7), we will calculate it 
using a Lorentzian distribution of A, only for numerical 
purpose, this is 



where 5 = 5i = [S^Jl) ^ = S,,. In Eq. (10) G^^,^. 



Pz(A) 



1 



TT 1A^ + Ab- 
using Eq. (8) in Eq. (7) we obtain 



(8) 



Gz{t + t/3) = e-*/*° e-^/^o. (9) 
with Iq as the characteristic decay time and tq = 3to- 

In Figure 3 we showed the measured results for the 
maximum values of the FID's, i.e. the pseudo-FID 

obtained by the values for an averaged amplitude close 
to t = 0, which are extracted by the reversion experiment 
in Figure 1 (c) using MREV8-like blocks of pulses as 
RS{i) in the isotropic 5CB (at T = 37° C). Besides, in 
Figure 3 the obtained data are fitted with an exponential 
decay function with characteristic time tq = 4.5ms. 
Comparing with the results of the same experiment 
in nematic LC's, whose measurements are shown in 
Figure 11 in QD-II, we can sec that for a reversion 
time of T ~ 2.5ms the signals in the nematic phase are 
vanishing but the decay due to inhomogeneity only can 
produce a reduction around of the 40% in such signals. 
Therefore, in the reversion experiment in nematic LC's 
the contribution to the decay of the signals from the 
field inhomogeneity can be neglected due to it is not an 
important effect affecting the dynamics. It is worth to 
note that it is possible to configure a new MREVS-like 
sequence introducing four tt pulses between the tt/2 
pulses in Figure 1 (b) with the aim of removing more 
efficiently the dynamics due to the field inhomogeneity. 
Nevertheless, such sequence would double of total time 
under pulses which adds non-ideal effects over the 
measurements like those who will be commented in in 
Section II of this supplementary material. Therefore, to 
avoid extra non-ideal effects and given that the decay 
produced by the field inhomogeneity has a longer time 
scale than the observed decays we used the standard 
MREV8 sequence. 

Finally, using Eq. (6) corresponding to the isotropic 
phase and the result for the FID signals under the re- 
version experiments on the nematic LC's, which was ob- 
tained in Section II. A in QD-II, it is ready to obtain the 
corresponding signals for an open quantum spin system 
in the nematic LC's under such reversion experiments, 
and affected by the field inhomogeneity, as 



X e-' Gc^,c; it) G (r) Gz{t + r/3), 

(10) 



and G, 



(rt) 



are the decoherence functions for the 



free-evolution dynamics and the reversion dynamics, 
respectively. It is important to note that these deco- 
herence functions depend on the eigenvalues C,^ and 
C( so imposing eigen-selectivity (see QD-II for details 
about eigen-selectivity). On the contrary, Gz does not 
depend on such eigenvalues thus it does not involve 
eigen-selectivity effects. 

We conclude in this section that the dynamics pro- 
duced by the field inhomogeneity has a time scale longer 
than the decoherence one and it does not involve eigen- 
selectivity. 



II. MISADJUSTMENTS IN THE DYNAMICS 
UNDER REVERSION 

In this section, we study with numerical simulations 
the effects of parameter misadjustments or errors in re- 
version experiments under the dynamics of a closed spin 
system. The aim of this study is to check for possible 
eigen-selectivity behaviours due to such misadjustment 
affecting only the spin dynamics. To do that, we use the 
eight-spins model of a FAA^e molecule, which was pre- 
sented in QD-I^. 

We simulated the dynamics under the reversion experi- 
ment of Figure 1 (c), where each block of reversion RS{i) 
is a sequence like the MREV8 one shown in Figure 1 (b) . 
To calculate the effect of the MREV8 sequence on the 
spin dynamics, first we define the dipolar spin Hamilto- 
nian of a nematic LC molecule as^ 



with the following 2nd order tensor 
T^2'o = 4 (3IzjIzk-^-I^) 



(11) 



V6 

1 



21zjlzk - 2 (I-f-jl-k + I-jl-l-k) 



(12) 



and the dipolar couplings defined as (in Hz units) 



_ fJ^oj'^h (l - 3 cos^ /3jk) 



87r2 



2r% 



(13) 



under the international system of units (SI), where the in- 
dices j, k run over all the proton sites within the molecule. 
In Eq. (13), 7 is the proton gyromagnctic ratio, /io is 
the vacuum permeability, Tjk 
tance between spins j and fc, 

lar angle of the vector r^t with respect to the system 
fixed to the molecule, Szz is the mean molecular order. 



is the internuclear dis- 
3jk stand for the po- 
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Ij = Ixj X + lyj y + Izj z is the vector of the i-th spin 
angular moment and I^i is the i-th spin angular moment 
operator, with a = {x, y, z, +, —}. Since all the spin op- 
erators in this section act on the spin Hilbert space of 
one molecule, for instance with the molecular index 1, 
for simplicity we will not explicit in such operators the 
upper index (si) pointing the corresponding molecular 
spin Hilbert space. 

Then, we are calculating a set of operators which result 
from applying pulses with an angle 9 and opposite phases 
on the dipolar Hamiltonian of one molecule, 'Hsi, such 
expressions will be useful in the calculation of the dynam- 
ics under reversion. The expressions are the following 



Ti'_i = ^ (IzjI-k + I-jIzk), (16b) 

T^2!V2 = ^I+jI+k, (16c) 

T^2% = il_jl_k. (16d) 

The Hamiltonian Hsi is called the secular part of the 

dipolar Hamiltonian and the Hamiltonians H^gl 

and "H^gl constitute the non-secular parts of the dipolar 
Hamiltonian. 



n 



(2) 
SI ' 



ill, {0) ^ lisi = ^ ^si 

- smecoseU'sl + d ^^-u'sl, 



3 cos2 e - 1 



2 



■Hsi 
3 sin^ 9 ^-^(2) 



2 2 



(14a) 



(14b) 



(14c) 



(14d) 



where the pulses of an angle 9 and direction a are defined 
as = e*'-", with I„ = J2i lai- In Eq. (14) we defined 
the Hamiltonians 



;i ^ttVqYI (^2'+! + T'2^i) , (15a) 



o7(i) _ 



Consider the case that the rf pulses in the MREV8 
pulse sequence of Figure 1 (b) are of an arbitrary 9 an- 
gle instead of 7r/2 pulses, due to experimental misadjust- 
ments. In such a case, and supposing the same angle 9 
for all the pulses, the evolution operator of the sequence 
in Figure 1 (b) will be 

Ri' U?^^T2/2)U(ri)Uj:,^(r2)U(ri)U?,^r2/2) R^ 
X R^ U?^^r2/2)U(ri)U^^^r2)U(ri)U^/(r2/2) Rt«, 



where 



U(ti) 



-iHsi Ti 



(17) 



(18) 



(19) 



with a = {x, y,x, y} and we applied the operator iden- 
tity 1 = RxRjc^ = Rl^Rx in tlie ends and in the middle 
of the sequence. 

When the times n and T2 are small enough, it is possi- 
ble to approximate the evolution operator (17) by joining 
together the exponent of all the operators (18) and (19). 
In such a case, and when the angle 9 is close to 7r/2, the 
contribution to the dynamics of the non-secular Hamilto- 
nians vanishes due to cancellation of the terms with T-L^gl, 

H^gl and H^gi, and the resulting evolution operator will 
be 



n^^l ^nVeYl ^ifc H'+i - T^2'-i) , (15b) 



ii'il ^ TT^e ^ Dj, (t^2^2 + T^2'_2) , (15c) 



with the following 2nd order tensors 



T2^+l — "2 (IzjI-l-k + I+jIzk) , 



(16a) 



ij'ksit2,h) Rl!\Jrt{t2)ij{h)K, (20) 
where we defined 

U(ii) = e"*'"^^*!, (21) 



Urt(t2) = e 



(22) 



with ti = 4ri, t2 = 4t2 and k = 2/(1 - 3cos^6l). It 
is worth to note that Eq. (20) has the form of the gen- 
eral evolution operator for reversion experiments anal- 
ysed in QD-II, except for the pulse operators at the ends 
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in Eq. (20). Nevertheless, the presence of such pulses 
does not alter the expectation values of the observables 
as calculated in QD-II; thus they can be neglected for 
simplicity in the study of decoherence effects in reversion 
experiments using the MREV8 sequence. If t2 = i^ti, 
we can see that \]'^jg{t2-ti) — 1- Therefore, for that 
time setting, the spin dynamics during ti will be per- 
fectly compensated (or will be very close to that) by the 
dynamics during the evolution time 1,2, and the initial 
state is recovered after the total time evolution. 
The pulse sequence of Figure 1 (c) consists of chains of 
blocks of MREV8 units. In this case, imder the same ap- 
proximation about the time parameters that yields the 
operator (20), we will obtain the same evolution opera- 
tor as in Eq. (20) for the experiment with n blocks RS 
of reversion, but with <i = 4nri and t2 = 4,nT2- 
In the case of an exact (and ideal) experimental setting 
of the 7r/2 pulses, the expressions in Eq. (14) simplify to 



nsi 



2 ^^1' 



(23a) 



■nsi 



-£7(2) 

"•SI ' 



(23b) 



and K = 2. 



To observe the effects over the dynamics under rever- 
sion due to the non-secular Hamiltonians and/or exper- 
imental misadjustments, we performed numerical sim- 
ulations of the experiment of Figure 1 (c), using the 
Hamiltonians (23) for a proton system model of a PAA^g 
molecule^, and including a parameter misadjustment of 
the reversion sequence, using t2 = {k + e) n (in this case 
K = 2). In the ease of small error in the setting of the 
tt/2 pulses, the value of k will be slightly different of 2. 
Besides, in Eq. (14) the factors 'sin ^? cos will be little 
different of zero and thus there will be a small contri- 
bution to the reversion dynamics coming from the non- 
secular Hamiltonians TL^gi and H^gi ■ However, if the 
time Ti is small enough (as well as T2) the contribution 
to the dynamics coming from the non-secular parts will 
be negligible and the; resultant effect will be a setting of n 
different of 2. In such cases, if we set T2 = 2t\, the error 
in the pulse angle will be equivalent to a misadjustment 
of T2 for a correct setting of the reversion sequence, due 
to K ^ 2; in this way we would be setting T2 = (k -|- e) ti 
(with K -|- e = 2). 

Since there are only two Hamiltonians in (23) with dif- 

ferent contributions of the non-secular Hamiltonian ngl , 
the sequence of operators given by Eq. (17) will actually 

be conformed by two equal sets of reversion evolution op- 
erators, i.e. there will be a duplication of such operators. 



Therefore, for the sake of simplicity we used in the sim- 
ulation of the reversion sequence, for each block RS, the 
following evolution operator 



USs(n) = Uj:j(K + e)ri]U(n)U?J(,^ + e)n]U(n), 



where 



U(ri) = e 



U^t [(k + e) n] = e-' ^si («+^) ^1 , 



(24) 



(25) 



(26) 



(27) 



with K = 2. 

The reversion evolution operator of the total n blocks has 
the following expression 

(") 

UH'sH(n) ^ U-^(ri) X ... X IJ^/siri) ^ U V^sin), 

(28) 

where the evolution operator is multiplied by itself n 
times. 

The evolution time of each reversion block RS is 



= {K + l + e)2Ti, 



(29) 



and the total evolution time has the discrete values of 
T = tiTc (n = 0,1,2,...). Thus, we can alternatively 
define the evolution operator (28) as 



Ufl'rw=u^'s(„)(n)' 



(30) 



where n = t/ [{k + 1 + e) 2ti]. 

The signals which we obtained from the simulations of 
the reversion experiments were the normalized expecta- 
tion values of ly and they were calculated as 



Iy(i,T)«;< 



ir,, {12} 



tr,, \ ly e 



-iHsit 



xU^'r(T)IyUj,«/'-(r)e^«-* . 



(31) 



We performed numerical simulations of Eq. (31) for dif- 
ferent values of ri and e with k. ~ 2. The results are 
shown in Figures 4-9. Figure 4 shows the behavior of 
Eq. (31) with the free-evolution time t and the rever- 
sion time r. Figure 5 shows the dependence with r of 
Eq. (31) for t = 0, which is useful to appreciate the rever- 
sion dynamics. Figure 6 shows the r dependent Fourier 
transforms of Eq. (31) in the time t; in this pictiue we 
can test the eigen-selection effects on the spectral lines 
of the signals. For such test, several cuts of the spectra 
for different normalized frequency lines as function of r 
are shown in Figure 7. 
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Figure 8 contanins the spectra of Figure 6 with the 
spectrum normalized for each value of t. In Figure 9 
several cuts of the spectra of Figure 8 for different 
normalized frequency lines are shown. With these 
normalized spectra it is possible to compare the decay 
rate in the time r for different spectral line. 

The evolution operator (30) can be approximated for 
small enough values of ri times as: 



U^'l'-(r) 



-i [2?isi + (K+e)(??gi+?i« J] nn 
-i [•Hsi-(1+€/k)«si] 2nTi _ JHsi ^^TTJ 



(32) 



where we used T-L^gi + "Hl^ = —2'Hsi/K, with k = 2, 
and 2nTi = r/ (k + 1 + e). Therefore, under this time 
approximation, using Eq. (32) in Eq. (31), we obtain 



(33) 



We can note that Eq. (33) is the normalized expression of 
(ly{t,Ty'^ shown in Section II. C in QD-II. The results 

showing the effects of the dynamics under the approxi- 
mation that yields (33) are those corresponding to the 
settings {n = 5/xs,e/K = 0.1} and {ri = 5/xs,e/«; = 0.3} 
in the Figures 4-9, where the effects of the non-secular 
Hamiltonians in the dynamics have been mitigated by 
adopting a small value of ti . We can see, using the trace 
property tr{AB} = tr{BA}, that for t = the signal ob- 
tained in Eq. (33) as a function of the reversion time r is 
the same as a single FID signal, but with a factor of time 
scaling of ■ Such kind of signals are shown in Fig- 

ure 5 for the mentioned setting, where it is observed that 
the behavior of the signals presents oscillations, which 
are typical of the FID's signals, where the smaller the ra- 
tio e/fc is, the lower the frequency of the oscillation will 
be. 

For general cases of values of ti and c/k, we can write 
the evolution operator (24) as 



UKs(n) = e 



(34) 



with 
basis 



Tc given in (29). 
{ I ai ) } common to 



We define the eigen- 
the evolution operator 



U^g(ri) and the Hamiltonian 'H|.^''^(ti), and we de- 
fine the eigenvalues j4j''^(ti) ~ (^ai |U^'^^(ti) | ai ) and 
ai''(Ti) = (ai |H^i''(ri)| ai ). Therefore the Hamil- 



tonian Tigi (ti) can be obtained writing its matrix 
representation in such eigenbasc, where Q!^''^(ri) = 
i lnAi''(ri)/Tc. Using Eq. (34) in Eq. (28), the evolu- 
tion operator (30) for the n blocks of reversion will be 



US5'"(T) 



in 



(35) 



with T = TiTc, and the expectation value (31) for general 
cases can be written as 



1 



trs. in] 



trs^ { ly e" 



lUsi t 



(36) 



The signals in Eq. (36) arc the normalized form of the 
expectation values (^Iy(f, r)^^ which were obtained in 
Section II. C in QD-II, while in Eq. (36) the dependence 
of 'hl}si with the parameters k, e and t\ has been 
explicitly written . 

It can be seen by comparing Eq. (32) and Eq. (35) that 
under the approximation of small values of n 



(37) 



and the expression of Eq. (36) is approximated by 
Eq. (33). 

In other calculated cases, if we set e = for dif- 
ferent values of ti, the signals (36) so obtained will 
show the effects on the reversion dynamics due exclu- 
sively to the non-secular Hamiltonians. The results 
which expose these cases correspond with the settings 
{ti = 10/is,e/K ~ 0} and {ri = 50/is,e/K = 0} in the 
Figures 4-9. In the extreme misadjustment setting with 
Ti = 50/is we can observe from Figure 5 a fast bell-shape 
decay of the signal, but such experimental setting can 
be improved drastically in the practice. On the other 
hand, such decay does not present an eigen-selection 
effect as can be observed from the comparison between 
the evolution in r of different spectral lines shown in 
Figure 6 and Figure 7 or the normalized spectra in 
Figure 8. Besides, in Figure 9 there is not a monotonous 
and eigen-selective decay for each spectral line with 
different frequency. 

Other results shown are for the settings 
{ri = 20jis,e/K — 0.1}, which constitutes an inter- 
mediate case between the ones with small values of ri 
and e ^ (expressed by Eq. (33)) and with large values 
of Ti and e = 0, and {ti = 5/is, e/k = 0}. The latter case 
would be close to the the best possible practical condi- 
tion and it is used as contrast with the other setting cases. 

The conclusion, extracted from the results of the sim- 
ulations, is that for the closed spin system there are no 
eigen-selection effects present in the signals produced by 
experimental misadjustments or contribution to the dy- 
namics coming from the non-secular Hamiltonians, as can 
be observed from the spectra evolution and the spectral 
lines evolution in the results shown in Figures 6-9. This 
conclusion reproduces the one extracted from the theo- 
retical analysis in Section II. C in QD-II. 
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^ H.H. Segnorile and R.C. Zamar, J. Chem. Phys. 135, 244509 FIG. 1: Experimental pulse sequences, a: Zeeman reversion. 
(2011) b: MREV8. c: Compound reversion sequence, where each 

reversion block RS{i) constitutes a single reversion sequence 
like the MREV8. 
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FIG. 2: Zeeman reversion experiment in nematic phase under 
resonance condition, a: 5CB (T = 27°C). b: PAA^e {T = 
115°C). 
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FIG. 3: Measured values of the maximum amplitude of the 
FID under a single MREV8 reversion, with continuous varia- 
tion of the setting times, for the 5CB at T = 37°C in isotropic 
phase under resonance condition (dots). The solid line is a 
fitting of the measured values with an exponential decay func- 
tion with a characteristic decay time of 4.5ms. 




FIG. 4: Simulated results for the FID's as a function of the 
free-evolution time t and the reversion time r, produced by 
reversion experiments with blocks of MREV8 sequence in a 
nematic PAAds molecule under resonance condition. The re- 
sults corresponding to different setting of the characteristic 
reversion time ri and the error time factor e/k are shown. 
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FIG. 5: Simulated results for the amplitude of the FID's for 
f = as a function of the reversion time r, produced by 
reversion experiments with blocks of MREV8 sequence in a 
nematic FAAds molecule under resonance condition. The re- 
sults corresponding to different setting of the characteristic 
reversion time t\ and the error time factor e/k are shown. 




FIG. 6: Simulated results for the Fourier transform in the 
free-evolution time t of the FID's as a function of the rever- 
sion time r, produced by reversion experiments with blocks 
of MREV8 sequence in a nematic PAA^e molecule under res- 
onance condition. The results corresponding to different set- 
ting of the characteristic reversion time ri and the error time 
factor e/k are shown. 
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FIG. 7: Simulated results for the evolution of different nor- 
malized spectral lines of the Fourier transform in the free- 
evolution time t of the FID's as a function of the reversion 
time T, produced by reversion experiments with blocks of 
MREV8 sequence in a nematic PAA^b molecule under res- 
onance condition. The results corresponding to different set- 
ting of the characteristic reversion time ri and the error time 
factor e/K are shown. 
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FIG. 8: Simulated results for the normalized Fourier trans- 
form in the free-evolution time t of the FID's as a function of 
the reversion time r, produced by reversion experiments with 
blocks of MREV8 sequence in a nematic PAAde molecule un- 
der resonance condition. The spectra ware normalized for 
each T value. The results corresponding to different setting of 
the characteristic reversion time ri and the error time factor 
e/K are shown. 
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FIG. 9: Simulated results for the evolution of different nor- 
malized spectral lines of the Fourier transform in the free- 
evolution time t of the FID's as a function of the reversion 
time T, produced by reversion experiments with blocks of 
MREV8 sequence in a nematic PAA^e molecule under res- 
onance condition. Where the spectra have been normalized 
for each r value. The results corresponding to different set- 
ting of the characteristic reversion time ri and the error time 
factor e/ft are shown. 



